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Experimental Study of Internal 
Gravity Waves Over a Slope 

by 

David A. Cacchione 

Submitted to the Department of Earth and Planetary Sciences 
in partial fulfillment of the requirement for the 
degree of Doctor of Philosophy. 



A series of laboratory experiments were conducted in a glass wave 
tank to investigate the propagation of internal gravity waves up a sloping 
bottom in a fluid with constant Brunt-Vaisala frequency. Measurements of 
the wave motion in the fluid interior were primarily taken with electrical 
conductivity probes; measurements in the boundary layer were made with 
dye streaks and neutrally buoyant particles. The results indicate that, 
outside of the breaking zone, the amplitude and horizontal wave number of 
the high-frequency waves increase linearly with decreasing depth; this is 
shown to agree with existing linear, inviscid solutions. A zone of breaking 
or runup is induced by these high-frequency waves well upslope. Shadow- 
graph observations show that, if the wave characteristics are coincident, 
or nearly so, with the bottom slope, the upslope propagation of the low- 
frequency waves causes a line of regularly spaced vortices to form along 
the slope. Subsequent mixing in the vortex cells creates thin horizontal 
laminae that are more homogeneous than the adjacent layers. These lam- 
inae slowly penetrate the fluid interior, creating a step -like vertical density 
structure. 

Available linear theoretical solutions for the velocity in the viscous 
boundary layer, determined to be valid for certain experimental conditions, 
are used to develop a criterion for incipient motion of bottom sediment in- 
duced by shoaling internal waves. The maximum sediment sizes that can be 
placed into motion, according to this criterion, are larger than certain 
mean sediment sizes on the continental margin off New England. This sug- 
gests that internal waves might induce initial sediment movement. Specula- 
tion about the geological effects of breaking and vortex instabilities is also 
given. These processes, not definitely measured in the field as yet, might 
also be conducive to sediment movement. 
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Assistant Professor of Geology 

Carl Wunsch 

Associate Professor of Oceanography 
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1. INTRODUCTION 



GENERAL 

The purpose of this research is threefold: 

(1) to provide experimental results that describe the changing 
character of single-frequency, small-amplitude internal gravity waves as 
these waves propagate over a sloping bottom; 

(2) to compare these results with recent theoretical solutions for 
those conditions which are theoretically tractable, and to examine the phy- 
sical characteristics of those conditions which are difficult to treat theo- 
retically; 

(3) to apply the theoretical and experimental results toward the 
development of a simple analytical model that prescribes the conditions for 
the initiation of bottom sediment motion induced by shoaling internal 
gravity waves . 

The motivation for this work was rooted in speculation - a starting 
point not uncommon in experimental research. The speculation was stimu- 
lated by the results of previous work that had indicated the possible signifi- 
cance of internal waves as geological agents and by recent theoretical 
solutions that predict the intensification of internal wave motion along 
sloping bottoms. 

It was decided early in this research that a clearer physical under- 
standing of the shoaling process for internal waves was needed - particuarly 
knowledge of the nature of the bottom boundary layer - before the two-phase 
problem could be considered. The effects of internal waves on sediment 
particles of various shapes, sizes, and densities were investigated during 
a preliminary set of experiments. It was concluded from these kinds of 
experiments that the reliability of the results and the range of test condi- 
tions were limited by the experimental difficulties, mainly the difficulty 
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of attaining large velocities in the small experimental system. The approach 
to this problem was changed to investigate the movement of sediment by in- 
ternal waves analytically, based on the results of a detailed laboratory study 
of shoaling internal waves over a smooth sloping bottom. Recent theoretical 
analyses had cast a framework for the experimental design (Wunsch, 1969; 
Keller and Mow, 1969); preliminary laboratory work had established the 
experimental techniques. Earlier experiments by other researchers had 
demonstrated the feasibility of some of the techniques (Mowbray and Rarity, 
1967; Martin, Simmons, and Wunsch, 1969; Gibson and Schwarz, 1963), 
and had shown several stimulating results (Thorpe, 1966). 

The subsequent work is divided into three sections: 

(1) the experimental design, instruments, and techniques are dis- 
cussed in Chapters 2 and 3; 

(2) the results of the experiments and their comparison with theo- 
retical solutions of Wunsch (1969) are presented in Chapter 4. The results 

* 

are separated into two basic types: experiments with high-frequency waves 
and experiments with low-frequency waves. The interior wave field and 
the boundary layer are discussed separately for each; 

(3) in Chapter 5 the balance of forces acting on a bed particle is 
used to derive a simple analytical model for the stability of the particle 
beneath a train of high-frequency, shoaling internal waves. The approach 
closely follows that used by Ippen and Eagleson (1955) and Eagleson and 
Dean (1959) in the development of an analytical model for sediment motion 



* 

The distinction between high and low frequency waves is discussed in 
more detail in Chapter 4. Basically the incident energy of the high fre- 
quency waves propagates upslope to the corner; i. e. , toward the intersec- 
tion of the slope and the free surface, with no reflection. By contrast, the 
incident energy of the low frequency waves theoretically is back reflected 
from the slope. 
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induced by surface waves. The applicability and limitations of the analysis 
of internal waves and sediment movement are discussed in Chapter 5, and 
some speculations are offered concerning the effects of breaking of internal 
waves and the action of boundary layer instabilities on bottom sediment in 
the ocean. 



BACKGROUND 

Internal Wave Theory and Experiment 

In two recent papers, Wunsch (1968, 1969) has developed normal- 
mode solutions for both standing and progressive internal waves in a wedge 
geometry. The solutions of interest here are for progressive waves pro- 
pagating into the wedge, i.e., over a linearly shoaling bottom. This theory 
assumes two-dimensional motion in a stably stratified, inviscid, Boussinesq 

2 £ 

fluid with constant Brunt- Vaisala frequency N (N = — — -• — ). The geometry 

P 0 dz 

and coordinate system that Wunsch considered is given in Figure 1. 

Wunsch’s solutions are discussed in more detail in Chapter 4; these solu- 
tions are presented for comparison with the experimental results. One of 
the principal features in the solutions is the suggestion of a strong intensi- 
fication of velocity components along the sloping boundary. Wunsch (1969) 
conducted a simple experiment in which internal waves were propagated 
over a linearly sloping, rigid bottom. From photographic determinations 
of the changes in wave length up the slope, he showed that the results are 
in agreement with his solutions. The wave length was approximately pro- 
portional to the depth. Magaard (1962) presented solutions for a similar 
problem of standing internal waves except for the density distribution and 
shape of the bottom profile. He considered a density distribution that can 
be expressed analytically as: 



p(z) = c 2 (2z + c x ) 



1/2 



( 1 - 1 ) 



where Cp c 2 are constants. The bottom profile that he chose for 
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computational facility and some realism was: 

1 / 2 “ 



h(x) = ^ 



1 



(2bx)‘ 



bx 



( 1 - 2 ) 



where c^ and b are constants. He developed a linearized characteristic 
equation for the velocity field and found that the equation was separable in 
the spatial coordinates only for the case of constant depth. He also showed 
that for a bottom slope y less than the slope of the characteristics c, a 
solution is possible which gives waves that have increasing amplitudes and 
decreasing wave lengths with depth. Magaard (1962) also found an appre- 
ciable intensification of the velocity components along the sloping bottom. 



Robinson (1970) has recently considered the effects of a corner 
(such as the point of transition between slope and flat bottom) upon an in- 
ternal wave train propagating toward the slope in an infinite medium (i.e., 
with free surface infinitely far from the horizontal bottom). Robinson 
provided solutions that consist of incident and reflected waves such that 
the reflected waves do not violate the radiation condition; i.e., they do 
not introduce energy from the far field (at x = + 00 ). 



Keller and Mow (1969) obtained an asymptotic solution to the 
problem of internal wave propagation in a horizontally stratified fluid of 
nonuniform depth by application of the principles of geometrical optics. 
Their general solutions apply to the linearized equations of motion, and 
the wave amplitude and fluid depth are assumed to vary very little over 
horizontal distances small compared to a wave length. They gave a partic- 
ular example in which N was assumed constant and the bottom was a small, 
linear slope. The solutions in this case show that the wave length is 
proportional to the depth, and the wave amplitude increases linearly with 
decreasing depth. Keller and Mow noted that their solutions also are 
in agreement with the experimental results given by Wunsch (1969). 

Hogg and Wunsch (1970; also personal communication) derived 
similar asymptotic solutions using a WKB approach with a two parameter 
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expansion scheme. Their first order, linear solutions for progressive 
internal waves over a small, linear slope are equivalent to those of Keller 
and Mow (1969); both solutions are compatible with the normal mode 
solutions of Wunsch (1969) for small slopes (y«c). Hogg and Wunsch ex- 
tended their analysis to higher-order terms in the stream function and 
developed expressions for the wave-induced radiation stress. They showed 
that a significant set-up of the isopycnals induced by internal waves is pos- 
sible and demonstrated theoretically the plausibility of "longslope” currents 
for a line of breaking internal waves oblique to the bottom contours. From 
the higher order terms Hogg and Wunsch (1970) were also able to demon- 
strate that in the absence of mixing the net Lagrangian motion of water 
particles in a two-dimensional stably stratified fluid is necessarily zero 
for internal wave motion in the proximity of a sloping, rigid boundary. 

Fofonoff (1969) also examined analytically the unsteady motion 
in a horizontally stratified, uniformly rotating ocean. He considered 
the seaward propagation of an internal baroclinic wave that was 
generated at the continental shelf edge by the diurnal tide. By exam- 
ining the characteristics of the solution, he presented a numerical 
example in which the values of density, latitude, and depth were those 
measured at site "D" (39° 20’N, 70°W), located about 37 miles south of 
the 200 meter contour. Fofonoff (1966) showed that north of this location 
the bottom profile of the continental slope and the path of the characteris- 
tic were approximately coincident. Again there is implication of increased 
bottom shear along the slope. 

Sandstrom (1966) derived the reflection properties of internal 
waves from a sloping, rigid boundary. Both he and Robinson (1970) have 
noted that the normal- mode solutions derived by Sandstrom from his ray 
theory solutions for progressive internal waves propagating over a slope 
violate the radiation condition (i.e., the complete solution requires an 
anomalous set of reflected waves that introduce energy from the far field). 
Sandstrom predicted that upon reflection the waves change their wave 
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number and amplitude; he found that the magnitude of the changes is 
dependent on a ratio equivalent to y /c. Sandstrom (1966) also conducted a 
set of experiments in which internal waves were generated by a flap-type 
wave maker and allowed to propagate over a planar, inclined bottom. His 
experimental results for a gently sloping bottom (y «c) showed qualitatively 
two features that were also found in the present study, and are discussed 
quantitively in Chapter 4: (1) amplification of internal waves traveling 
from "deep" into "shallow" water, and (2) intensification of motion near 
the slope. Sandstrom (1966) noted, however, that unfortunately the partic- 
ular experiment for which these features were observed was complicated 
by the presence of a surface mode that was also generated by the wave 
maker. Although his apparatus was apparently capable of generating the 
first internal mode, there are no experimental results in his work which 
indicate that he did this. Sandstrom discussed solutions for the case y >c 
(steep bottom slope), but did not present experimental results for this 
case. He concluded from the linear theory that in the case y ^ c, the motion 
close to the bottom becomes very large. 

Longuet- Higgins (1968) extended the discussion of reflection prop- 
erties of internal waves from rigid surfaces to account for various types 
of bottom roughnesses. He has shown that the important parameter is not 
the ratio of the scale of bottom roughness to wavelength but instead the 
ratio of the scale of roughness to the thickness of the oscillatory boundary 
layer. He pointed out that this suggests that small-scale irregularities 
can severly affect the transmittance and reflection conditions of the bottom. 

In a very thorough theoretical and experimental treatment of finite- 
amplitude effects on interfacial and internal waves in a horizontal channel, 
Thorpe (1968) also reported some qualitative observations of progressive 
internal waves on a slope. In his experimental arrangement, a linearly 
sloping, smooth bottom was mainly employed as a wave absorber. How- 
ever, during the course of the experimental runs for internal waves in a 
linearly stratified salt solution he noted several interesting features over 
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the slope. These observations are summarized here as a preface to the 
results shown later in this study: 

(1) the wave length decreased upslope; 

(2) a large runup was apparent along the upper slope region, and the 
larger amplitude waves produced a larger runup; 

(3) the free surface remained undisturbed; 

(4) a "rotor" was observed near the free surface on one occasion; 

(5) breaking ("overturning") was observed well upslope on one 
occasion; 

(6) no reflections could be detected. 

Thorpe (1968) showed three photographs of the waves over the 
slope. His method of flow visualization, involving alternating dyed and 
undyed layers, did not permit detailed observations of the breaking or 
mixing, but his qualitative observations are quite instructive and essentially 
agree with the results shown later. 

It should also be mentioned here that excellent experimental studies 
of instabilities and turbulence in a stably stratified fluid, including breaking 
of interfacial waves, have been carried out by Pao (1968). 

Field Observations 

Two kinds of field observations are pertinent to this study: (1) 
measurements that suggest shoreward propagation of internal waves over 
the continental margin, and (2) observations that indicate, directly or in- 
directly, the motion of bottom sediment due to passage of internal waves. 

Measurements of the first kind have been presented by Lee (1961), 
Gaul (1961) and Ufford (1947), among others. Their investigations presented 
time-series measurements of temperature fluctuations taken simultaneously 
at three stations oriented in various triangular arrays in shallow water 
(depth <200 feet). The wide geographical distribution of these studies illus- 
trates the ubiquity of internal motions on the continental shelves. Common 
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internal- wave periods were on the order of 5 to 20 minutes; wave propaga- 
tion was shoreward with speeds of 0.3 to 1.2 knots. Lafond (1962) has 
provided similar values. Longer period internal oscillations are described 
by Lee (1961), Boston (1964), and Summers and Emery (1963). The last 
study illustrated the refraction of internal waves of semidiurnal period as 
they propagated shoreward over the continental slope and shelf off Southern 
California. The estimated wave speed was 7 knots in deep water and slightly 
less than 1 knot over the shelf. 

Isolated asymmetrical internal temperature disturbances resembling 
the idealized solitary wave profile of surface waves (Ippen, 1966) have also 
been detected by Lee (1961), Gaul (1961), and Cairns (1967). The common 
period of this sort of disturbance seems to be approximately 10 minutes, 
and the propagation direction is again shoreward. 

Curves of spectral kinetic energy density versus frequency are 
available for velocity measurements of internal motion taken with vertical 
arrays of current meters at site "D” (Fofonoff, 1968). An example of these 
spectra is shown in Chapter 5. 

Direct observations which reliably link the two motions of internal 
waves and sediment movement are limited to a few isolated instances. 

This lack of evidence might suggest both the difficulty and the insufficiency 
of this kind of field measurement. Ideally, one would like to measure simul- 
taneously the internal wave velocity structure and the movement of bottom 
sediment beneath the wave motion. Either measurement alone is a challenge 
in the field. Lafond (1965), using underwater television to observe the sedi- 
ment motion and a vertical array of temperature sensors to measure the 
internal wave motion, found that the movement of sediment ripples in shal- 
low water (depth about 60 m) was correlated with large fluctuations in the 
temperature structure. He interpreted this as evidence for sediment trans- 
port by internal waves. 

Indirect field data that connect the two motions are similarly limited. 
Revelle (1939) inferred the connection in an attempt to explain an observed 
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pattern of sediment distribution. Munk (1941) calculated the theoretical 
standing-wave characteristics in the Gulf of California and found that the 
nodal-antinodal separation coincided favorably with Revelle's (1939) ob- 
served distribution of sediment sizes. That is, the smaller sizes were 
noticeably lacking beneath the nodes, where larger instantaneous velocities 
are expected. However, in a subsequent detailed study of the sediment dis- 
tribution in the Gulf of California Van Andel showed that the sediment pat- 
tern described earlier by Revelle was not apparent in the later data. 

Emery (1956; 1960) suggested that deep, standing internal waves in certain 
California basins stir the bottom sediments and might produce the observed 
local concentrations of coarse sediments at the basin sills. 

More recent investigations of the large submarine sand waves (av- 
erage height =* 8 meters) near the edge of the continental shelf southwest of 
Great Britain (water depths 180 meters) by Stride and Cartwright (1958), 
Cartwright (1959), Stride and Tucker (1960), and Carruthers (1963) indicate 
the possibility of interaction between internal waves and bottom sediment. 
Cartwright (1959) formulated a theory which suggests that these sand waves 
are maintained by a stationary internal- wave system of tidal period that 
exists during the presence of the local seasonal thermocline. Measure- 
ments of oscillations of the scattering layer (Stride and Tucker, 1960) and 
towed thermistor records (Carruthers, 1963) support this theory. Car- 
ruthers (1963) also measured the currents close to the crest of the sand 
waves at one-hour intervals over a seventeen- hour period. A plot of cur- 
rent vectors indicated that the predominant flow direction was normal to 
the crests, with a maximum current of 0.5 knot in this direction. Inter- 
estingly, the crests of the sand waves are parallel and adjacent to the 
shelf edge. 

Many investigators have studied the observed sediment distribution 
on the continental shelf and slope off the East Coast of the United States. 
Uchupi (1963), Emery (1966), and Schlee (unpublished manuscript) have 
summarized much of this work. Sand waves and local topographic highs 
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are abundant on the continental shelf. The type and size distribution of 
sediment varies both across and along the bottom contours. The seaward 
variations are generally sharper. For instance, Uchupi (1963) delineated 
four broad zones parallel to the shelf edge from Cape Cod to Hudson Can- 
yon. He noted that the most seaward zone (about 70 kilometers wide) con- 
sists of silty sand, sandy silt, and silt (water depth 60 to 135 meters). 
Uchupi (1963) claimed that from Hudson Canyon to Cape Hatteras the outer- 
most zone is composed of fine sand of approximately uniform grain size. 
The surface of the continental slope in these regions consists mainly of 
silt and clay; however, sand, silty sand, and gravel line the floors of the 
submarine canyons that cut the slope. 

Wave-Induced Bottom Sediment Motion 

The writer could find no published results on wave-induced bottom 
sediment motion in a stratified fluid. However, there are many published 
studies of sediment motion beneath surface waves in a homogeneous fluid. 
Most of these studies are empirical or semi-empirical. The obvious lack 
of theoretical treatment probably stems from the complex form of the gov- 
erning equations of motion, which are not well established. 

One appealing approach to a laboratory and analytical study of the 
problem was undertaken by Ippen and Eagleson (1955) and later extended by 
Eagleson, et al (1958), and Eagleson and Dean (1959). A statistical treat- 
ment of laboratory data that involved the motion of discrete spherical sedi- 
ment particles due to shoaling surface waves supplemented their theoretical 
presentation. Many of the complicating secondary factors, such as (1) 
nonuniform bed roughness, (2) shape variations of the particles, (3) mu- 
tual particle interactions, (4) nonuniform incident wave trains, (5) local 
channelization and fluidization of the bed, and (6) bed permeability were 
not included in their work. They were able to derive a force-balance 
equation for the individual particles. Ippen and Eagleson (1955) described 
two possible types of particle motion: incipient and equilibrium. Incipient 
sediment motion was defined by Eagleson et al (1958) as "an instantaneous 
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condition reached when the resultant of all active forces on the particle 
intersects the line connecting the bed particle contact points. " They also 
defined established sediment motion as "an oscillatory or quasi- os dilatory 
condition of motion reached when for some portion of each wave cycle the 
sum of the instantaneous active forces is greater than that value necessary 
to initiate motion." The major difficulties in working with these definitions 
are: (1) the point of application of some of the individual forces is usually 
not known; (2) the various coefficients, such as those for drag, lift, and 
virtual mass, are not well defined for unsteady motion; (3) usually the 
effective fluid velocities acting on the bed particle must be estimated from 
an approximate theoretical equation, since the exact velocity profile near 
the rough bed in unsteady motion is not known. Despite these difficulties, 
Eagleson and Dean (1959) were able to arrive at a semi-empirical relation- 
ship for particle size, beach slope, and local wave characteristics . Eagle- 
son et al (1958) compared the analytical results with field data from eight 
stations on the Atlantic coast of the United States. As these authors put it, 
the agreement was "fortuitously" remarkable. 

The concept of an equilibrium diameter for given wave conditions 
was first proposed by Cornaglia (Miller and Zeigler, 1964) and has since 
been labeled the null-point theory. His arguments, as well as those of 
Ippen and Eagleson (1955) and Miller and Zeigler (1964), are basically 
simple. In an area of shoaling waves, every bottom sediment particle is 
affected by the existing sea state. For any given sea state and bottom 
slope, there will be one size class of particles in a state of equilibrium 
motion. That is, this size class will move back and forth along the bottom 
as successive troughs and crests pass, but will oscillate about some mean 
position. In general, particles larger than the null-point size will move 
offshore and those smaller will move onshore. Miller and Zeigler (1964) 
noted that in the field the null point should be regarded statistically, where- 
as in the laboratory it is convenient to regard it as single-valued. Their 
field study also indicated the limitations and the modified applicability of 
this theory in nature. However, at the same time their work showed that 
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the theory gives a satisfactory prediction of the gross sediment distribution 
of a nearshore sand bottom. 

Other laboratory investigations have produced relationships between 
the initiation of motion of a given particle size and the associated wave con- 
ditions (Manohar, 1956; Goddet, 1960; Vincent, 1958). These investigations 
were all conducted in a channel with a level bottom. Vincent (1958) observed 
that the fluid velocity just outside the boundary layer was almost constant 
at the onset of movement of equal- diameter grains. Abou-Seida (1965) 
noted that field measurements of radioactive sand tracers conducted by 
Sato et al (1963) compared favorably with the laboratory relationships de- 
duced by Manohar (1956) and Goddet (1960). 

A second approach to this problem is relatively new. Arguing from 
basic physical principles and certain assumptions, Bagnold (1963, 1966) 
derived a relationship between bed load transport and the dissipation of 
energy in the fluid motion. His argument rests on the assumption that a 
distribution of dispersive and tangential shear stresses exists in a bed of 
granular material that is subjected to shear deformation. Several sum- 
maries of the topic of wave-induced sediment motion are available (Inman, 
1963; Abou-Seida, 1965; Raudkivi, 1967). Such presentations of incipient 
sediment motion usually include a discussion of the various criteria developed 
for open channel flow. In particular, the work of Shields (1936), White 
(1940), and Kalinske (1947) is often cited. 
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2. EXPERIMENTAL SYSTEM 



APPARATUS 

Wa ve Tank and Mechanical Components 

The rectangular wave tank illustrated in Figure 1 has an open top 
and 3/8-inch glass side walls and bottom. It is 4. 9 m long, 20. 5 cm wide, 
and 38. 4 cm deep and consists of two sections of similar dimensions joined 
end to end by an O-ring seal contained between plastic inserts attached to 
the sections with epoxy cement. A thin layer of clear silicone sealant (Dow 
Corning 701 Building Sealant) was applied along the inside seam to form a 
smooth transition between the two sections. All other joints were also 
formed with silicone sealant. The ends of the channel are acrylic plastic 
plates sealed by Neoprene gaskets. A 14-foot length of aluminum angle 
attached to five vertical aluminum angle supports behind the tank and 3/4 of 
an inch above the side wall served as an instrument mount. 

A filling technique similar to that described by others (Fortuin, 

1960; Oster, 1965; Jaffee, 1968) was designed to establish the desired sal- 
inity gradients in the wave channel. Figure 1 shows a diagram of the filling 
system. Owing to the corrosive properties of salt water, the insides of the 
two 55 -gallon containers were coated with epoxy paint, and all seams were 
caulked with silicone sealant. Wherever the liquid came into contact with 
bare metal, only stainless steel or brass was used. Flexible Tygon tubing 
is used throughout the filling system wherever tubing was required. 

The filling rate was controlled by a brass needle valve. Each drum 
was equipped with a Tygon sight tube for visual determination of the level 
in each barrel. A motor driven stainless stirring rod with two brass stir- 
ring vanes was mounted above the fresh water barrel. To prevent solid- 
body rotation of the water, the stirring rod lay off the axis of the barrel. 

The bottom slopes, one of which is shown in position in Figure 1, 
were cut from Plexiglas sheets. Each slope is beveled at one end to form 



tank 6. Tygon filling tube 
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a smooth junction with the bottom of the tank. Four slope angles were used, 

7, 15, 30, and 45 degrees. During the experiments the slopes were supported 
by two parallel Plexiglas struts. These struts were positioned at the proper 
angle along the inner walls of the tank. Four Lexan plastic rods with threaded 
ends and plastic nuts were used to separate the struts and to maintain their 
position by forcing them against the channel walls. Neoprene strips were 
cemented to the outer surfaces of each strut to form a continuous and flexible 
seal with the side walls of the channel. All materials forming the beach 
slopes and supports were noncorrodible. 

Wave Generating Apparatus 

The wave generator had to meet five primary requirements: 

(1) smooth and sufficiently powerful motion; 

(2) desired frequency range and control; 

(3) variable amplitude; 

(4) compact size with components stably mounted on a plate above 
one end of the channel; 

(5) long running with low vibration level. 

In an effort to meet these requirements the apparatus was modified several 
times. Availability of various components such as the variable speed drive 
often aided the choice of parts. The end product was an outgrowth of the 
usual experimenter’s dilemna: getting optimum performance at reasonable 
cost. 

Pow er unit. Figure 2 is a diagram of the power drive for the wave maker. 

The entire drive mechanism is rigidly supported above one end of the chan- 
nel byfour 1-1/2-inch pipe legs with foot flanges. The base is a 1/4-inch 
steel plate. The assembly consists of a Graham variable- speed drive 
coupled to a 1/4-hp motor. The output shaft is joined to a fixed-ratio gear 
reducer (9.5:1) via a flexible coupling. This second speed reducer is con- 
nected to a Scotch yoke device that converts the rotary shaft movement into 
rectilinear sinusoidal motion of a vertical shaft. 
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Wave maker . Figure 3 shows the wave maker in actual location. Several 
additional considerations were involved in the design of this device. First, 
a unit was needed that could be inserted into position easily without appre- 
ciably disturbing the stratification after the tank had been filled. Secondly, 
the driving motion had to generate, to a good approximation, first- mode 
internal waves while minimizing mixing at the wave maker. In addition, 
the water motion in the section behind the wave maker had to be prevented 
from interacting with the fluid in the working section of the channel. 

Basically, the device consists of a plastic flapper plate free to pivot 
about its midpoint. Hinges anchor this plate to a supporting frame. A thin 
flexible rubber gasket covers the entire face of the plate and frame, thus 
preventing liquid from flowing between the two sections . This gasket is 
held in position by thin plastic covers cut to the same shape and dimemensions 
as the frame and flapper plate. The flapper plate is driven by a stainless- 
steel push rod connected to the vertical shaft of the Scotch yoke by a universal 
coupling that allows for minor misalignments between shaft and flapper plate. 

A closed-cell Neoprene strip was cemented to the sides and bottom 
of the frame. This afforded an adequate seal between the wave flapper 
assembly and the glass tank. A test showed that the wave maker could suc- 
cessfully maintain a head difference of at least 6 cm of water between the 
working section of the channel and the unused portion. 

DATA ACQUISITION 

This section describes the instruments used to record the experi- 
mental data. A more detailed discussion of the quality of the measurements, 
including calibration curves, error, and reliability estimates, is presented 
in the next chapter and in Appendix A. 

Detection of Wave Motion 

The changing properties of the internal waves as they propagated 
upslope were measured in the Eulerian sense with conductivity probes. 
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Fig. 3. 



Wave maker in actual location, 
depth is 38.2 cm. 
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The basic feature of the water medium was its strong salinity stratification. 
Also, over the range of salinity values used during this study, salinity is 
approximately linearly proportional to conductivity. The temperature vari- 
ations in the water were found to be small (usually less than 0. 5°C over the 
total depth) except in the upper two centimeters, so that below this layer the 
total change of conductivity at a point is here a valid estimate of the total 
change of salinity: 

s = s 0 + L - V 

(2-1) 

Da _ 1. DS 
Dt y Dt 

where 

Tj in %o/(ohm - cm)" 1 
S in %o 

• u - 1 - 1 

a in ohm cm 

Because the wave-induced fluctuations of conductivity are measured relative 
to some basic state of rest, only relative changes in conductivity at a point 
are needed to decode the waves. 

The conductivity sensor is a spherical platinum electrode that is the 
active arm of an A. C. Wheatstone bridge. Design and construction of the 
probe was inspired by Gibson and Schwartz (1963). Figure 4 shows a probe 
assembly, including electrode, glass housing, and connecting cable. A de- 
tailed description of its construction and an analysis of its performance are 
presented in Appendix A. 

Temperature 

Temperature of the room air and the water in the filling barrels 
was measured with a mercury-in-glass thermometer which can be read to 
0.05°C. Vertical and horizontal temperature profiles of the stratified 
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water column were taken with a glass-coated thermistor bead that is the 
active arm in a D.C. Wheatstone bridge. 

Specific Gravity 

Specific gravity measurements of the water in the filling barrels 
were made with a lead- weighted glass bulb hydrometer. The hydrometer 
can be read to 0.0005 for any specific gravity in the range from 1.000 to 
1. 070. Its accuracy is about 0. 2 percent of the maximum value in the range. 

The basic specific gravity gradient was determined from measure- 
ments of the vertical distribution of index of refraction. Figure 5 shows 
the refractometer and mini-siphon that were used to obtain index of re- 
fraction at various depths in the stratified water. These measurements 
were converted to specific gravity values from calibration curves. The 
siphon consists of a 0. 5 mm diameter stainless steel tube supported verti- 
cally by a brass bushing and two locking screws. The bushing is press- 
fitted into an aluminum bar. Centering pins allow repeatable positioning of 
the aluminum bar athwart the channel. Markings that are evenly spaced one 
centimeter apart are inscribed on the steel siphon for vertical positioning. 
Small-diameter Tygon tubing fits over one end of the siphon tube. A brass 
pin that can be inserted into the open end of the Tygon tubing is used to stop 
the flow. Flow rate through the siphon at mid-depth (about 16 cm of water) 
is about 0. 2 m £ /sec. 

Wave-Mak e r Motio n 

A rotary potentiometer in a D.C. voltage -divider circuit transduces 
the oscillating, rectilinear motion of the Scotch yoke into a fluctuating D. C. 
voltage signal. A precision gear meshes with a linear rack that is cut into 
the forward part of the upper vertical shaft of the Scotch yoke (Figure 2). 

The spur gear is fastened to the shaft of the rotary potentiometer. Vertical 
movement of the Scotch yoke causes a change in the voltage output of the 
sensing circuit. The sensitivity is approximately 0. 1 volt/cm. This tech- 
nique establishes a convenient and accurate frequency reference and time 
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Fig . 5 . 



Ref ractometer and mini-siphon. 
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base for the conductivity probe measurements. The accuracy of the circuit 
as a frequency reference is primarily dependent upon the linearity of the 
potentiometer. In this case, the linearity of the rotary pot is 0. 1 percent 
(Giannini Control Corporation). Amplitude calibration of the signal also 
provides a measure of the forcing amplitude of the wave flapper. 

The period of the motion is sensed by a microswitch mounted on 
the lower part of the vertical support plate for the Scotch yoke (Figure 3). 
This switch is closed once per cycle by a protrusion on the rotating arm, 
and provides a convenient trigger for a period counter. 

Photographic Measurement 

Photographic measurements of the wave motion were obtained by 
recording on film the behavior of neutrally buoyant particles and various 
dyes. Expandable polystyrene particles (Sinclair-Koppers Company) were 
chosen because their density was adjustable to the proper range by a tech- 
nique described by Bohlen (1969). In addition, their milky white color pro- 
vides an excellent contrast with the surrounding water when used with a 
black background and overhead lighting. After processing to achieve the 
desired density range, Bohlen found that their sizes varied over a consider- 
able range. Sieving produced a more uniform size distribution, with median 
size about 0.2 mm. The range of specific gravity for the lot was 0.990 to 
1. 040, with an obvious skew toward the heavier limit. The mean specific 
gravity was estimated at 1.020. These values were adequate for neutral 
buoyancy in the range of fluid densities used during this study. The particles 
were used to study orbital velocities and modal structure. 

Potassium permanganate (KMnO^) was used to make the vertical dye 
streaks to determine near-bottom velocities and bottom shear stress. 

When dropped into a water column, crystals of this substance leave a dis- 
tinct dark- red dye trace from surface to bottom. The vertical trace was 
then tracked photographically as it was deformed by wave motion. 
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Small quantities (one to two cc) of common liquid food coloring or a 
highly nondiffusive blue dye (Blue Dextran 2000) were injected into the fill- 
ing tube by means of a 16-gauge hypodermic needle and syringe through a 
self- sealing rubber membrane at infrequent intervals. The membrane was 
fitted over one end of a glass T-joint in the Tygon filling line (Figure 1). 

These injections resulted in thin horizontal dye layers at various levels in 
the stratified water column. This enabled qualitative measurements of the 
internal waves by observation of the deformation of the dye layers by the 
propagating waves. The layers were typically 3 to 5 mm thick for the food 
dyes and only about 1 to 2 mm thick for the Blue Dextran 2000. The tech- 
nique also aided in observing flow conditions at the entrance to the main 
channel during filling. Excessive mixing at this point was indicated by rapid 
diffusion of the dye, as might be expected in an area of strong turbulence. 

In such cases, filling rate was decreased. Finally, the dye layers added 
some aesthetic value to the final stratification. 

Shadowgraph images were used to display instabilities in the boundary 
layer and the violent breaking of waves in the upslope turbulent zone. The 
shadowgraph system is shown diagrammatically in Figure 6. The two-way 
diffusive screen is sold by Edmund Scientific Company under the trade name 
Lenscreen. Its excellent dispersive properties and good transmittance pro- 
duced clear and evenly lighted images. The merits and techniques of shadow- 
graph flow visualization are discussed by several authors (Barnes, 1954; 
Goldstein, 1965; Edgerton, 1958). A Bolex 16 mm camera was used to take 
motion pictures of the shadowgraph images. 

The shadowgraph lighting was provided by a tungsten filament lamp 
(Beck Company) that contains a dial rheostat for light-intensity adjustments. 

It also has a circular portal lens that can be moved toward or away from 
the light source. This lens and an adjustable diaphragm are located in a 
collimating tube through which the light rays pass. A 500-watt combination 
flood- spot lamp provided illumination for the dye streak photographs. 
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Fig. 6. Schematic of shadowgraph system. 
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Data Recording and Display 

Figure 7 shows some of the instruments used in this study. The 
oscilloscope at the left of the picture has provisions for dual-channel opera- 
tion (Tektronix Type 564 Storage Oscilloscope), and was used mainly to 
display the output of each conductivity bridge when null was set (Appendix A). 
The instrument at the right is a Brush dual-pen strip-chart recorder. The 
two charts operate on the same time base but with independent voltage scales. 
This recorder provided a real-time display of the transduced conductivity 
fluctuations during the internal-wave runs. Thermistor temperature mea- 
surements were also recorded in this way. 

The heart of the data recording system, a Wang model 2300 digital 
readout and recording system (Wang Laboratories), is at the center of the 
photograph in Figure 7. The fastest available sampling rate (1 cycle per 
second) was chosen for all phases of this investigation. The system was 
programmed such that each tape record was in card-image format (< 80 
characters per record) and represented one second of real time. Included 
in each record were digital clock time to the nearest second and six sequen- 
tial channels of data. The data characters contained five measurement 
digits, decimal point location, and sign; each data point was recorded to the 
nearest millivolt. The recording format simplified use of the data in the 
digital computation routines on the IBM 360/67 system and in listing the 
raw data in decimal format on the IBM 1401. 

Several cameras were used to photograph the neutrally buoyant par- 
ticles and dye lines. General photographic work was done with a Nikon F 
single-lens 35 mm reflex camera. Another Nikon F camera with an auto- 
matic film advance attachment was used to photograph deformations of the 
dye streaks, and, in some cases, the motions of neutrally buoyant particles. 
The film advance was controlled by a synchronous timing device that provides 
continuous one-second trigger pulses to the camera at an accuracy of ±0. 1 
percent of the timing interval. A Calumet view camera with fast Polaroid 
film was normally used to photograph the neutrally buoyant particles, which 
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Fig. 7. Data acquisition system 

Left: Tektronix type 564 stroage oscilloscope 

Center: Wang model 2300 digital data system 

Right: Brush dual-pen strip chart recorder 
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were illuminated by a flashing stroboscope. The strobe flashes were triggered 
externally by a signal generator (Wave-Tek). 

EXPERIMENTAL PROCEDURE 

Although the procedure for each experiment varied depending upon the 
purpose, the following basic procedure was employed throughout the prepara- 
tion and execution of all experiments. Several days (usually three) prior to 
the experiment, the two 5 5 -gallon barrels were filled to a predetermined level 
with hot tap water to which 7 to 10 teaspoons of Alconox laboratory detergent 
were added. Liberation of dissolved air accompanied the subsequent cooling 
to room temperature. The detergent inhibited formation of air bubbles along 
the walls of the drums. The water was also stirred manually as often as pos- 
sible. On the next day, a predetermined amount of salt was added to one drum. 
A submersible pump was used to mix the salt into solution and to minimize 
entrainment of air at the free surface. After the wave tank had been cleaned, 
the wave flapper assembled and tested, and the slope supports positioned at 
the desired angle, the siphon connecting the two barrels (Figure 2) was primed, 
and the stirrer in the fresh water barrel was started. The water temperature 
in each barrel and the room temperature were then recorded. Readings of 
specific gravity and index of refraction were taken in both barrels. During 
filling, the water level was recorded as a function of time, and dye was in- 
jected into the filling tube at selected intervals. Room and water temperatures 
were monitored at one-hour intervals. After filling was completed, about 10 
hours later, the final temperature and specific gravity of the water in each 
barrel were measured. The slope was then inserted very slowly and smoothly 
by sliding it down the supporting struts until its beveled edge lay flush with 
the tank bottom. Very little mixing was observed, except at the edge, as the 
slope moved through the stratified fluid. The coherence of the dyed layers 
near the slope attested to the small degree of disturbance. The wave flapper, 
whose bottom had been positioned at mid-depth during the filling, was then 
inserted through the remainder of the fluid. Its driving rod had been fastened 
to the vertical shaft of the Scotch yoke prior to filling. This procedure also 
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caused very little noticeable mixing. The water depth was recorded, and 
the top of the tank was covered with a thin plastic sheet (Saran Wrap) to 
minimize evaporation losses at the free surface. (Prior experience showed 
that these evaporative losses could be quite substantial, 1-2 mm of water 
per day. ) In addition, the conductivity probes were cleaned and platinized 
as necessary (Appendix A). Depth markings were checked and refurbished 
on all probes. 

On the next day, vertical profiles of temperature and index of re- 
fraction were taken over the flat bottom and at other stations along the slope. 
The depth was recorded, and the presence of an interface near the surface 
was investigated with the shadowgraph apparatus. The plastic covering was 
removed from the slope section of the tank, and the pertinent photographic 
apparatus was made ready. 

Prior to wave generation, the conductivity probes were given static 
calibration by raising and lowering them vertically by known increments 
while recording their output at each level; this allowed a comparison with 
previous calibration values and provided an updated sensitivity coefficient 
for each probe. Normally, the next group of measurements was taken 
with the probes all at the same depth over the flat-bottom region. Adjacent 
probes were separated by equal distances usually 5 to 10 cm, with the probe 
farthest from the slope at a distance of at least 200 cm from the wave maker. 
Runs were then conducted in which the forcing amplitude and frequency of 
the wave maker were varied to include most of the values to be investigated. 
These data provided a set of reference values with which the later measure- 
ments over the slope were compared. The procedure used to take the dye- 
streak measurements is summarized below. 

(1) Prior to the initiation of wave motion, a grid marked in milli- 
meters was gently lowered to the slope in the middle and parallel to the side 
walls of the tank, photographed in place, and gently removed. The position 
of the camera relative to the front wall of the tank was carefully measured; 
this position could be repeated to within 1 percent. 
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(2) After the disturbances introduced by the grid removal had sub- 
sided (usually 10 - 15 minutes), and several minutes after the wave-maker 
had been turned on, several dye pellets were dropped at various locations 
over the slope (normally four) and at one position over the flat bottom. The 
pellets were dropped through a small glass tube at the previous position of 
the grid. By choosing dye pellets that were small (diameter <0.3 mm) and 
nearly spherical, irregular motion during descent was minimized. During 
the rest of the experiment, photographic and conductivity measurements 
were taken at various locations and times. Each experiment usually lasted 
from two to four days, depending on the nature of the particular experiment 
and the physical endurance of the experimenter. 
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3. DATA ANALYSIS 

Several kinds of experimental data were obtained; these included 
(1) temperature, (2) mean density gradient, (3) density perturbations, 

(4) boundary layer velocities and bed shear stress, and (5) miscellaneous 
data including mode number, net particle drift, and measurements of boun- 
dary layer instabilities. The instruments and some of the procedures that 
were used to obtain the various kinds of data have already been described; 
this section presents a summary of the data reduction techniques and the 
methods of analysis. 

TEMPERATURE 

Vertical temperature profiles at several horizontal locations were 
made at least once before and after each experiment, and usually several 
times in between. Typically the maximum temperature variations were 
relatively small (<0.5°C) over the fluid interior; more irregular changes 
often occurred in the near-surface layer. This layer will be discussed in 
more detail in the next section. Figure 8 is a representative temperature- 
depth profile. 

MEAN DENSITY GRADIENT 

A typical density gradient is shown in Figure 9. Measurements of 
index of refraction were made at about 2 -cm intervals in the vertical and 
at closer spacing near the free surface. Index of refraction was converted 
directly to specific gravity by applying a calibration constant that included 
a correction for the local temperature. Several times during the course of 
the experiments the specific gravity values obtained in this way were checked 
by carefully weighing a known volume and determining density in g/cc di- 
rectly. The maximum difference in the two methods was always found to be 

3 

less than 1 part in 10 (or about 0. 1 percent). 



0 

2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

22 

24 

26 

28 

30 

"32 

£ 

Lq. 



32 



Z , cm 




1 1 1 1 1 1 — 

.5 23.0 

T emperature , 



t r 



° C 



. Typical temperature variation 
Data from exp. group 4. 



23.5 



th depth. 



0- 

2- 

4 

6 - 

8 

io- 

12 

14 - 

16 - 

18 - 

2Q 

22- 

24 - 

26 - 

28 - 

30 

3 2 



33 




ig. 9. Typical specific gravity variation with depth. 
Data from exp. group 2; computed vertical 
distribution of Brunt Vaisala frequency is 
also shown. 
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The linearity of the curve in Figure 9 is quite good over the fluid 
interior but nonlinear tails are evident near the bottom and the free surface. 
Evaporative cooling at the free surface, when the plastic wrap was removed 
from the slope section of the tank during the conductivity runs, tended to 
increase the density of the water there, countering the effect of lower den- 
sity due to low salt concentrations. The overall result was a more homo- 
geneous density layer near the free surface. Weak convective motions were 
observed in this layer in the absence of wave motion by dropping potassium 
permanganate crystals through the water and tracking the motions of the dye 
traces in this layer. The amount of evaporation was reduced significantly 
(to < 1 mm per day) when the plastic cover was used. 

The variation of salt concentration with depth as a function of time 
for a similar experimental arrangement has been derived theoretically 
from the one- dimensional diffusion equation for a constant coefficient of 
diffusion and conditions of no flux of salt across the boundaries (Mowbray 
and Rarity, 1967). The theoretical concentration-depth curves were shown 
to have a linear interior and nonlinear tails that thickened in time. The 
point of maximum curvature located between the nonlinear and linear por- 
tions of the curves might explain the sharp shadowgraph image that was 
usually observed a small distance below the free surface and was inter- 
preted as a change in the density gradient (i.e., ap Q /a z + 0 at this point). 

An example of the vertical distribution of Brunt-Vaisala frequency N is 
shown in Figure 9. The effects of the nonlinear tails cause N to decrease 
sharply in the small layers. Photographic measurements of the modal 
structure and of the particle displacements indicate that the basic wave 
motion was not measureably affected by the small, nonlinear portions of 
the density distribution. 

DENSITY PERTURBATIONS 

Perturbations in the density field induced by internal wave motion 
were sensed with conductivity probes. Voltage variations relative to some 
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mean or null value and proportional to the conductivity changes were recorded 
digitally on magnetic tape. The quality of the data was checked by using a 
real-time display on a strip-chart recorder. 

Normally, measurements from all four probes were taken at least 
one minute before and after the wave motion to establish an instrument drift 
correction. Each of six data channels was sampled sequentially once per 
second (At = 1 sec); the channel allocations were four channels for probe 
data, one for the wave -maker signal, and one for a reference voltage. 

Samples of the raw data signals from one probe and wave-maker are shown 
in Figure 10. The length of each data record depended on the particular ex- 
perimental run. Record durations of least 50 wave periods, often up to 100, 
were obtained for quasi-steady conditions as determined from the strip- 
chart output. 

The data records on the magnetic tape were then processed on the 
IBM 360/65 as follows: 

(1) average and trend were computed, listed, and removed; 

(2) mean square of the data values was computed and listed; 

(3) the periodogram was computed for 2N data points using a Fast 

Fourier Transform algorithm (IBM Manual H20-0205-3) to obtain squares 

of the real Fourier coefficients a, , b, such that: 

k’ k 



N-i 

v = l/2a + ^ l\ cos (^) + b k sin(^)} + ^a N (-l) j 
J k=l 

j = o, 1 , , 2N-1 (3_1 ' 

Normally, the starting point for each computation of the periodogram was 

selected when it was determined that the data were quasi- steady (transients 

and other irregular motions were no longer evident on the strip chart). 

Output listings and selected plots of the squares of the periodogram esti- 
2 2 

mates (in cm ) and of the phase at each frequency f (in cps) 
were obtained. 
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Fig. 10. Raw data plotted for wave maker and conductivity 
probe output. Mean not removed. (Data from 
Exp. Group 5) 
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m = k = 0, 1, N-l 

A sample periodogram computed from the data record represented in Fig- 
ure 10 is shown in Figure 11. A calibration constant has been applied so 

2 2 

that the ordinate is cm instead of (volts) . The strong peak at the input 
frequency and the sharp fall-off at adjacent frequencies are typical of these 
results. The cut-off frequency was lower than the Brunt- Vaisala frequency 
for all periodogram computations (cut-off frequency = 0.5 cps). 

2 * 

Estimates of the periodogram a Q were computed for probe data 

taken over the flat bottom region. These estimates were used as reference 

2 

values for subsequent estimates a for data taken over the slope at similar 

m 2 2 
frequencies. The results are shown in a later chapter in the form a^ /a 

to indicate the change in wave amplitude from flat bottom to slope regions. 

This normalization also conveniently removes certain calibration constants. 



VELOCITIES AND BOTTOM SHEAR STRESS 

Time- sequenced photographs of the deformation of potassium per- 
manganate streaks were used to obtain one-second averages of the Lagran- 
gian velocities of these streaks at several levels parallel to the bottom in 
and near the bottom boundary layer. Reliable measurements of streak dis- 
placements were obtained at one-second intervals over at least one full 
wave period, often down to distances of 0.05 cm above the bottom. Tracks 



* 

The subscript "o" signifies a value pertinent to the flat-bottom section of 
the tank, not to be confused with the general Fourier coefficient a in 
Equation (3-1). 
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of neutrally buoyant floats near the bottom generally showed very little, if 
any, motion normal to the bottom, except in the zone of breaking. The dye 
traces usually remained coherent to the slope, but if instabilities and tur- 
bulence were present the streak traces usually developed a wavy structure 
or dissipated rapidly near the bottom. 

The following summarizes the methods used to obtain velocity and 
bed shear stress data from the streaks. 

(1) Photographic negatives which contained a series of one- second 
positions of the streaks for a particular input wave amplitude and frequency , 
and at each station, were enlarged and transferred onto a sheet of tracing 
paper by using a microfilm reader. A single sheet was produced containing 
the traces of the successive positions of the dye streak at each station for 
one- second sampling intervals. Successive displacements of neutrally 
buoyant particles were recorded simultaneously on the same sheet. Figure 
12 is an example of a worksheet that includes the dye streak and particle 
positions. The numbers represent successive seconds in time relative to 

t = 1 . The streaks and particles were also color-coded; this does not show 
in Figure 12. 

(2) The coordinates of each streak relative to its position at the slope 

were punched onto IBM cards by a special purpose X-Y digitizer. The streak 

positions were usually sampled at 0.02 mm increments normal to the bottom 

* 

in the boundary layer and at larger increments (0.05 mm to 0. 10 mm) above 
the boundary layer. Each neutrally buoyant particle position was digitized. 

(3) Successive displacements of the streaks (or particles) at a partic- 
ular level above the bottom gave approximately one- second averages of the 

Lagrangian velocities u. , at that level; 

t 



* 

The thickness of the boundary layer is tentatively defined as the distance 
of the inflection point in the streak trace to the bottom. 
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t = 1, 2, , T 

j = 1 , 2 , , L 

T is the number of streaks (usually equal to the nearest integer 
greater than the wave period). 

L is the number of vertical levels. 



j = 1 represents the lowest level for which the data of any particular 
set are considered reliable. 



As stated above, the spacing between the various vertical levels 
from j = 1 was constant in the boundary layer. 

(4) Estimates of bottom shear stress (one-second averages) were 
computed from the velocity values at j = 1 ; 



t = 1, 2, , T as before. 

A 77 is the distance between the bed and the level j = 1 measured normal to 
the bottom. 



MISCELLANEOUS 

To ascertain that the wave maker was in fact generating primarily 
the first mode, several photographs of the orbital motion of neutrally buoy- 
ant particles over an entire vertical section were taken at various positions 
along the tank. Figure 13a is a typical section over the flat bottom. The 
predominantly vertical motion near the center of the water column and the 
horizontal pattern near the free surface and at the bottom are characteris- 
tic of internal wave motion of the first mode. The photograph was taken 
with Polaroid film that was exposed by flashing a strobe through a thin slit 
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Fig. 13 a. Time sequence of neutrally buoyant 

particle positions over horizontal bottom. 
Modal structure approximately equal to the 
first internal mode (n = 1) . 

Pointer depth: 16.0 cm. 

Bright horizontal section on pointer: 1.0 cm. 
Brunt-Vaisala period: 6.5 sec. 

Wave period: 12.0 sec. 

(Strobe: 1 flash/sec; Polaroid 4x5 film, 

ASA 3000) 
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from above the tank at 2 cps. The illuminating flash was a rectangular 
beam about 1. 5 cm wide and approximately centered between the side walls 
of the tank. Similar photographs were taken at various positions over the 
slope. Figure 13b shows a close-up view of successive displacements of 
neutrally buoyant beads near the bottom of a relatively low slope ( a & 15 
degrees). Note the amplification in the particle displacements nearer the 
bottom. 

Shadowgraph images were used to visualize the boundary layer in- 
stabilities and breaking of the high frequency waves. Both of these phenomena 
are discussed in Chapter 4; several photographs of the shadowgraph images 
are shown in that chapter. 

The net drift patterns of neutrally buoyant particles were determined 
by two methods: 

(1) The positions of selected particles were marked on the glass side 
walls of the tank over long durations (up to 30 minutes). 

(2) The positions of the particles were traced from 16 mm movies 
of the shadowgraph images that were projected onto a screen of large grid 
paper. By counting the frames between position marks, the approximate 
time that had elapsed was estimated and a net velocity computed. 
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Fig. 13 b. Time sequence of neutrally buoyant 

particle positions over slope (a=30 deg) . 
Note intensification of motion near bottom. 
Brunt-Vaisala period: 6.0 sec. 

Wave period: 9.6 sec. 

Input wave amplitude: 0.2 cm. 

(Strobe: 1 flash/sec; Tri-x 35mm film.) 
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4. EXPERIMENTAL RESULTS 
GENERAL DISCUSSION 

This section describes the organization of the results and presents 
an overview of the hydrodynamics. Subsequent sections discuss the results 
in more detail and compare them with the predictions of recent theories. 

The results are separated into two basic categories: (1) one that describes 
the interior wave field (wave number and wave amplitude); (2) another that 
describes the activity at the boundaries (wave breaking and boundary layers). 

It is useful to define two distinct hydrodynamic regimes on the basis 
of the reflection conditions for internal gravity waves from a smooth, rigid 
slope (Phillips, 1966; Sandstrom, 1966). This subdivision is based on the 
ratio of bottom slope y to the slope of the input wave characteristics c. 
Figure 14 illustrates the classification and shows that regime n can be fur- 
ther subdivided into two cases. Phillips (1966), among others, has shown 
the slope of the characteristics to be a function only of the stability fre- 
quency N and the frequency of the input waves co, provided that N is con- 
stant with depth: 

c = 1/(N 2 /o> 2 -1) 1/2 (4-1) 

Since N and y were essentially constant during each experiment, the ratio 
y/c could be altered by changing In this way, both hydrodynamic regimes 
y /c< 1 andy/c>l were investigated during several experiments. 

Tables Bl and B2 in Appendix B summarize the experimental param- 
eters for the various measurements that were made with the conductivity 
probes. Table Bl is divided into experimental groups that each have a simi- 
lar Brunt-Vaisala frequency and bottom slope. This table is a reference for 
Table B2, which lists the various input wave frequencies that were generated 
during the experiments. The experimental conditions for each frequency can 
be found by matching the second number in the first column of Table B2 with 
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Reflection 



Regime 


Case 


Designation 


y /C 


Properties 


I 


A 


subcritical 


< 1 


transmissive; all 
energy propagates 
into corner 


n 


B 


critical 


= 1 


marginally trans- 
missive; energy 
propagates parallel 
to bottom slope 




C 


supercritical 


> 1 


nontransmissive; 
energy is back- 
reflected 




Fig. 14. Classification of internal waves based on reflection 
from a smooth, rigid slope. Slope of wave charac- 
teristics = c; bottom slope = y. 
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the number of the experimental group in Table Bl. Those frequencies 
which belong to hydrodynamic regime II are identified by a dash in column 7. 
An asterisk is located above those frequencies for which more than one ex- 
perimental run were made. 

The emphasis on two hydrodynamic regimes is supported by the re- 
sults. The subcritical case was examined extensively, since the results for 
this case are pertinent to the later analysis of sediment motion and since 
this case is important for certain oceanic conditions (Chapter 5). The 
data for this case suggest a further classification into two wave types 
based on the frequency of the input waves. This distinction is more justified 
for the measurements that were taken over the low slope angle (15 degrees) 
than those measurements taken over steeper slopes. The experimental re- 
sults indicate that the high-frequency input waves decrease in wave length 
and increase in wave amplitude as they propagate upslope until at some lo- 
cation near the corner and below the free surface they dissipate by breaking. 
The low-frequency waves in this case (y<c) initially steepen like the high- 
frequency waves, but the effects are less pronounced. The steepening 
diminishes upslope as these waves form an internal surge that produces a 
large runup with little breaking. The water motion near the bottom is am- 
plified for both wave types. The bottom shear stress is increased over the 
slope and is greatest in and near the zones of breaking and runup. In addi- 
tion, there is considerable vertical motion near the bottom in the breakers. 

The experimental results for case B and case C (y/c>l) indicate 
that the interior wave motion is very complicated. The conductivity-probe 
measurements show that during most experimental runs the wave amplitude 
over the slope was a maximum at or near a depth z = -cx. (The geometry 
is defined in Figure 15.) In the critical case, a line of small, regularly 
spaced vortices forms along the slope. The vortices have diameters on the 
order of 1 cm and axes that are normal to the side walls of the tank. They 
oscillate along the slope with the wave motion, collapsing and reforming 
over a half-cycle. The distance from the vortex core to the slope is about 
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Fig. 15. Geometry and coordinate system. 
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7 to 10 mm. Several wave cycles after the formation of the line of vortices, 
thin horizontal layers appear in the fluid. Each layer is associated with a 
vortex near the slope and represents a modification to the density structure 
over the slope. Wave breaking near the corner accompanies the vortex 
activity. These features are discussed in more detail later in the chapter. 

THEORETICAL DISCUSSION 

The experimental results presented in this chapter are compared to 
the theoretical solutions of Wunsch (1969). Preliminary measurements in- 
dicated that the exact linear solutions for the interior wave field fitted the 
experimental data. The experimental conditions and the assumptions im- 
plicit in the theory were compatible except for the boundaries. Wunsch (1969) 
solved the problem for a wedge geometry bounded above by a rigid horizon- 
tal top and below by a rigid bottom slope that extends indefinitely far from 
the corner (the term corner refers to x = z = 0 in Figure 15). The rigid 
bottom slope used in the experiments terminated a finite distance from the 
corner at the transition to a horizontal bottom. The possible effects of 
this discrepancy between theoretical model and experimental conditions 
are mentioned in the discussion of results for cases B and C and in the 
boundary layer section for case A. Although a free surface was permitted 
in the experiments, the observed vertical motions of this surface were 
small. In the sense that w « 0 at z = 0, the use of a free surface was 
compatible with the rigid top in the theoretical model. 

Fora viscous Boussinesq fluid of mean density p Q (z) the two-dimen- 
sional perturbation equations (after Wunsch, 1969) are: 

u t = -P x /p 0 + ^ v2u ( 4-2 ) 

w t = -P Z /P 0 - g p7p 0 + W 2 w (4-3) 



(4-4) 
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p' t + W P Q = 0 (4-5) 

— . z 

where p' is the perturbation density, and a bar denotes partial differen- 
tiation. 



Interior Wave Field 



The interior equations are given by the inviscid forms of Equations 
(4-2) through (4-5) above. The pertinent linear solution (Wunsch, 1969) 
for small-amplitude internal waves of frequency a> traveling toward the 
corner is composed of two progressive waves: 



ip = 



A 



e ~iqfn(cx - z) 



e ~iq£n(cx + 




(4-6) 



2n7 1 
q “ £n A 




(4-7) 



where n is the mode number. This solution is singular at x = z = 0 for 

case A (y <c). The phases of the two waves in Equation (4-6) are arith- 
metically averaged to give an average phase function U: 



U = q/2 [£ n(cx - z) + £n(cx + z)] 



(4-8) 



It follows that an average local horizontal wave number is defined by k : 



k 

w 



9U/9X 




(4-9) 



The corresponding wave number k^ that was determined from the conduc- 
tivity probe measurements is A0 m /Ax (Chapter 3). 

• 1 

For a time dependence of e w , the incompressibility 
condition (Equation (4-5)) yields: 



p = -i/w 8p 0 /az w 



(4-10) 



(the prime on p has been dropped for convenience). 
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The vertical velocity w is obtained from the stream function in Equation (4-6), 
w = a i ///ax: 



w = Aiqc 



1 -iq£nM 
M e + 



1_ -iq£nN 



-icot 

e 



(4-11) 



M = cx - z , N=cx+z 

Applying Equation (4-11) into Equation (4-10) and taking the time average of 
the product of p and its complex conjugate p*, we obtain an expression for 
the square of the amplitude variation with coordinate position: 




1^ = | A 2 ( 3£ ) 2 (8Po/82) 2 



X 




2 

NM 



cos {qCnkS} 



M 



(4-12) 



The comparative experimental value is the square of the spectral amplitude, 
2 

a , as determined from the periodogram of the conductivity-probe measure- 

m 2 2 
ments. The method of comparison between a and p was discussed in 

m w 

Chapter 3. 



It is useful to point out here that fory<<c, Wunsch's solutions are 
the asymptotic solutions for internal waves propagating over a small slope 
(also given by Keller and Mow, 1969). The dependence on local depth h 



for local values of horizontal wave number k and amplitude a is 



k = 
x 



n7rc 



a = 



a h 
o o 



(4-13) 



where the subscript "o" signifies "deep" water values (i. e., constant values 

of wave amplitude a and depth h ). Hogg and Wunsch (1970) have shown 

o o 

that the nonlinear terms are significant when e becomes order one, 



a h 

, o o 
e = ak = — r — cnji 
x h 



(4-14) 



That is, the amplitudes of the next higher order term in the velocity ex- 

2 

pansion grows as e . This suggests that the linear solutions are no longer 
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valid over the small slope at some depth h given by e of order one in 
Equation (4-14). 

For the case y>c, Wunsch (1969) showed that the apparent singu- 
larity at z = -cx in Equation (4-6) could be removed. Wunsch rewrote q 
* 

in the form 



q = r + is 



(4-15) 



and showed that the expression for the stream function for waves propagat- 
ing upslope could be rewritten 



i/i = A 



(s - ir)fnM 
e 



e (s - ir) i nN 



-icot 

e 



(4-16) 



where 

_ 2n7rgn|A / 1 
[(fn|A'|) 2 + i, 2 ] 



-2n7r 

[ ( C n | A |) 2 + t r 2 ] 



A! = S-Z.2 

C + y 



1/A 



and n is a negative integer. (Note that £nN is complex for cx + z < 0.) 



5jc 

The notation here does not conform to that of Wunsch (1969) in order to 
avoid later confusion. He wrote the complex form of q as q = rj + is, 
and gave an analogous expression for Equation (4-16) in which the real part 
of the exponential group was written in a power of s(i. e., s in his notation). 
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In order to have the velocities and their first derivatives continuous, 
Wunsch derived the requirement that s > 2 or 



-n > 



(£ n 1 A' |)' 



7T 



(4-17) 



-n represents the lowest, nonsingular inviscid mode permissible for 
these solutions. The vertical velocity is obtained from Equation (4-16): 



w = Ac(s - ir) 



M 



e (s - ir)PnM 



\ 

N 



e (s - ir)£nN 



- i cot 
e 



(4-18) 



2 

An expression for p w analogous with that in Equation (4-12) can be obtained 

in principle from Equation (4-10) and Equation (4-18) and compared to the 

2 

measured amplitude a^ for specific coordinate positions by the method 
described in Chapter 3. However, it is not simple in this case to obtain a 
physically meaningful phase function that is analogous to Equation (4-8) 
owing to the complicated wave forms in Equation (4-18). The amplitude 
dependence on the spatial coordinates indicates that these are not simple 
plane waves. By writing 



w = f(x, z) e 



ig(x, z, t) 



(4-19) 



where f and g are complex, a phase function U can be defined as 

U = tan' 1 $$ (4-20) 

where I(w) and R(w) are the imaginary and real parts of w, respectively. 
This form of U will be used in the later discussion of theoretical and ex- 
perimental results for cases B and C. 



Boundary Layer 

The full set of perturbation Equations (4-2) through (4-5) were used 
by Wunsch (1969) to obtain linear boundary- layer solutions for case A 
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(y < c). The solutions match those for the inviscid interior just outside 
the boundary layer. The stream function can be written: 

^T “ + ^b.0. (4-21) 

where ^ and ^ are the interior and boundary layer contributions, re- 
spectively. In terms of the coordinates J, rf shown in Figure 15, the 
expression for ^ ^ given by Wunsch (1969) can be rewritten 



Vi. 



F e -iqCn(Bf) _ Q e -iqfn(DD 



x e -(i + 1 ) rf/5 e -iwt 



(4-22) 



where 



P = 



i - 1 2 

\w 



2 coi/ \ 
- N^sin^oy 



v V2 



6 = (i - 1) P 
B 



= c cos a + sin a, 



F = 



c sm a ~ cos a 
B ’ 



D = c cos a - sin a 



G = 



c sin a + cos a 

D 



The component of velocity parallel to the slope (in the 4 direction) is 
given by 



Urp — Uj + U b ^ — “d^rj^/dTj 



(4-23) 



or from Equation (4-6), after it has been expressed in £, tj coordinates, 
and Equation (4-22), 
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U T = V (?) 



l _ e "(i + !) rj/6 



-icot 



V(?) 



_ 



Fe -iqtn(Br) _ G e ~iq£n(DJ) 



(4-24) 



The behavior of the velocity component parallel to the bottom with 
variations in J and ? is obvious: (1) u^ grows to infinitely large values 
as ? -* 0; in fact, at if = j - 0, the solution is singular, as might be ex- 
pected since the interior solution is singular at the origin for y < c; (2) 
for if = 0, £ * 0, the solution reduces to that of the interior field; and (3) 
as "f -*> + oo, the boundary-layer part (u^ ^ ) becomes infinitesimally small 
and u T — vu. The solutions for w^ can be derived in principle from 
w T = a ip/dj, but are not presented here since the measured values of w^ 
were very small in the boundary layer (this is discussed later in this 
chapter). 

For cases B and C (y > c), Wunsch (1969) found that the boundary- 

layer equation for ip^ was nonseparable in the spatial coordinates. In 

fact, Wunsch showed that when the slope is near critical the scale of the 

-1/3 -1/2 

boundary-layer thickness is proportional to R ' , in contrast to R for 
y <c, where R is a local wedge Reynolds number and h is the local depth: 

R = Wc?/ V 



INTERIOR WAVE FIELD - RESULTS 
Subcritical Case 

Table B3 (Appendix B) summarizes the periodogram amplitudes 
and wave numbers that were determined from the conductivity-probe mea- 
surements as described in Chapter 3. The notation in the column headings 
is discussed in the following sections. The experimentally determined 
values are organized into columns by input wave amplitude. 
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Wave number . As was described in Chapter 3, an experimental wave 
number was computed from the measured phase differences between 
adjacent probes (k = A0 m /Ax). Columns 9, 10, and 11 of Table A3 con- 
tain the appropriate values of divided by q, i. e., ^ A0 m /Ax (q is de- 
fined in Equation (4-7)). This parameter is convenient since ^k ffl corres- 
ponds to the theoretical value ^ given by Equation (4-9): 



1 k = l au = X 
q w - q 0 X - x 2 . (z/c) 2 



(4-25) 



The right-hand side of this equation is evaluated for the appropriate values 
of c, and the average is computed for the two probe positions. The value 
listed in column 7 is 

i k i 



2 / , v2 

x - (z/c) 



1 9 . 



(4-26) 



in which the subscripts 1 and 2 indicate that the function is evaluated at the 
coordinate positions for each probe and then the average is taken. 

The theoretical curves in Figures 16 and 17 are straight lines with 

slope of 1. The data in Figure 16 are taken from the values for relatively 

high-frequency waves in Table B3. The upper right portion of the diagram 

corresponds to the corner region of the slope. Figure 17 is a similar plot 

for the data of the low-frequency waves. Both plots show that the wave 

* 

number decreases with increasing distance from the corner. The agree- 
ment of the results with the theoretical curve is best for waves over the 
lower section of the slope, well away from the corner. Toward the corner, 



The behavior of the abscissa, 



x - (z/c y 

constant z. (This is obviously true for x >> z/c.) 



, is approximately like 1/x for 
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Fig. 16. 

Ups lope variation of wave number, y < c 
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Fig. 17. 

Upslope variation of wave number, y < c. 
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the data points in Figure 16 depart significantly from the theoretical curve. 

It is also evident that the value farthest to the right for each experimental 
group in Figure 16 is displaced toward higher wave number than that pre- 
dicted by theory. This trend is illustrated in Figure 18, which includes 
data for the high-frequency waves (co = 0. 094 cps) in Group 1 (Table B3). 

The upslope increase in wave number is similar for both input waves 
(a Q = 0.2 cm and a Q = 0. 1 cm). However, the greater increase in wave 
number for the larger waves (a Q = 0. 2 cm) at the most cornerward position 
(i. e., closer to x = 0) suggests that the upslope variation in wave number 
might depend on amplitude. The experimentally determined value of e at 
this most cornerward position is approximately 0.15; this might indicate 
that nonlinear effects are beginning to have a significant effect on the motion. 

In Figure 17 there is good agreement between theory and measure- 
ment for the low-frequency waves, although most of the data values tend 
to be lower than the theoretical curve in the right-hand side of the diagram 
(near the corner). This suggests that the wave length of these lower fre- 
quency waves decreases more slowly with decreasing distance to the corner 
than theory predicts. 

Another type of comparison between theoretical and experimental 

values of wave number is shown in Figures 19 through 21. The data are 

taken from columns 12 through 14 in Table B3 (Appendix B). Experimental 

wave numbers k m and k Q were determined for each input wave from probe 

measurements taken over the slope and over the flat bottom, respectively. 

The ratio k^/k^ is the ordinate in Figures 19 through 21 (k^ = A0 m /Ax). 

A comparative theoretical ratio is used as the abscissa in these figures; 

i. e., k w /k WQ where k w is the theoretical wave number computed from 

Equation (4-9) for coordinates established by the probe locations over the 

slope, and k is similarly computed for coordinates at the entrance to 
w o 

the slope region (k = A0 w /Ax). (The depth z used in the computation of 

k is equal to the probe depth during the determination of k and k Q .) 
o 
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Fig. 18. 

Upslope variation of wave number, y < c. 
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Fig. 19. 

Measured versus computed values of wave number, y < c. 
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Fig. 20. 

Measured versus computed values of wave number, y < c. 
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Fig. 21. 



Measured versus computed values of wave number , y < c 
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Figure 19 and 20 also show better agreement between experiment 
and theory for the smaller wave numbers over the lower slope sections. 
Farther upslope the data for the higher frequency waves in Figure 19 tend 
toward larger values of wave number relative to the theoretical values; by 
contrast, the lower frequency waves in Figure 20 generally tend toward 
smaller wave number. The effect of input wave amplitude on wave number 
is indicated in Figures 21 and 22, for high and low frequency waves, re- 
spectively. The values for the larger amplitude waves in Figure 21 appear 
to deviate from theory more significantly than the smaller waves. Signifi- 
cant departures from theoretical estimates are also shown for the larger 
input waves of low frequency (Figure 22). In this case, however, the 
measured wave numbers are smaller than the predicted values at locations 
farthest to the right in Figure 22. This behavior corresponds to the pre- 
vious observations (Figure 17) that the wave lengths of the lower frequency 
waves apparently decrease less than theory predicts at positions well upslope. 

Wave amplitude. The results of the conductivity-probe measurements for 
case A, concerning the effects of shoaling on wave amplitude, are sum- 
marized in Table B4 (Appendix B). The experimental group numbers cor- 
respond to those in Table Bl. Experimental determinations of the period- 

2 2 

ogram amplitudes a m and a Q for each input wave were made from probe 

measurements over the slope and over the flat bottom, respectively (for 

more detail see Chapter 3). The value a was determined for each of the 

2 2 0 

four probes so that each ratio /a Q represented data taken with a par- 
ticular probe. This procedure was convenient since it removed the 
necessity of introducing calibration constants relating one probe to another. 
That is, the measured amplitude ratios could be compared to one another 
without applying relative gain factors. This method of normalization also 

simplified the detection of measured changes in wave amplitude relative to 

2 

the input wave amplitude a Q . The comparative theoretical values p w and 
p ^ were computed from Equation (4-12) for coordinates of the probe loca- 
tions over the slope and at the entrance to the slope region, respectively. 
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Fig. 22. 



Measured versus computed values of wave number, y < c. 
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2 2 

The ratio p^ /p Q was equally convenient, since its use circumvented a 
conversion from density p to the units of a and permitted direct com- 
parison with the experimental ratio a ^/ a Q- 

The results shown in Figures 23 and 24, for the high-frequency and 
low-frequency waves, respectively, suggest that the agreement between 
theory and experiment is best for the lower values of ( a m / a 0 ) > i. e., for 
the wave motion on the lower and middle sections of the slope. This is in 
agreement with the results for wave number. The solid lines in these fig- 
ures were drawn for reference and represent points where measurement 
and theory coincide. The deviations of the points from the straight line in 
Figure 23 are significant for (p /p ) >20, but no consistent pattern is 
apparent in the variations. The values for lower frequency waves show 
earlier departures from theoretical values than do the higher frequency 
waves. For example, the measured and theoretical values shown for low- 
frequency waves in Figure 24 are significantly different for (p /p ) > 10. 
There is strong indication in the trend of the points that the actual ampli- 
tudes are much lower than predicted on the middle and upper sections of 
the slope. At these slope sections it was also found that the measured 

values of steepness a k (or e ) are considerably less than those esti- 

m m nr 

mated from theory. 

Measured amplitude growth is further compared with theoretical 

estimates in Figures 25 through 30. In these figures the variation of wave 

amplitude is examined as a function of horizontal distance x fi to the corner 

at various fixed depths z . The data are again taken from Table B4 (Appen- 
* c 

dix B). In each of these figures the theoretical and experimental ratios of 



* 

The actual wave amplitude in centimeters at a point along the slope can 
easily be found from the graphs. If the correct wave amplitude (in centi- 
meters) for the flat bottogi region is multiplied by the square root of the 
measured value ((a /a )), the measured estimate of the actual amplitude 
is obtained. m 0 
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Fig. 23. 



Measured versus computed values of wave 
amplitude, y < c. 



68 




Fig. 24. 

Measured versus computed values 
of wave amplitude, y < c. 
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Fig. 25. 

Upslope variation of wave amplitude, y < c. 
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Fig. 26. 



Upslope variation of wave amplitude, y < c. 
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Fig. 27. 

Upslope variation of wave amplitude, y < c. 
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Fig. 28. 



Upslope variation of wave amplitude, y < c. 
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Fig. 29. 



Upslope variation of wave amplitude, y < c. 
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Ups lope variation of wave amplitude, y < c. 
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2 2 

amplitude variations (p w /p 0 ) and ( a m / a 0 ) > respectively, are plotted 
along the ordinate; the ratio of the horizontal distance between the probe 
and the corner x fi to the total horizontal projection of the slope x q is plotted 
along the abscissa. The depth of the probe measurements is constant for 
any one of the figures, but as shown in Appendix B, Figure B4, it generally 
varies for each experimental group. In Figures 25, 26, 27 and 30 the re- 
sults pertain to the slope angle of 15 degrees, for which the distinction be- 
tween high and low frequency waves appears most applicable. Two effects 
of the method used to display the results must be mentioned here: 

(1) the deviations between the theoretical and measured estimates 
of amplitude are exaggerated in the graphs, since squares are plotted as 
the ordinate; 

(2) the peaking of the waves is a function of the vertical position z 

c 

of the probes as well as distance from the corner This latter effect is 
subtle; since the internal waves are of the first mode, the amplitude should 
be largest near the center of the water column at any position over the slope 

jjc 

if the modal structure remains intact. Consequently, the amplitude is 
expected to vary along a horizontal line (constant depth) due to the modal 
structure as well as to the amplification during shoaling. 

To the right of the peak values (i. e., away from the corner) in 
Figures 25 and 26 (high-frequency waves), the variations are essentially 
parallel, although the experimental values are consistently higher than the 
predicted. The maximum measured values are slightly higher than pre- 
dicted, with an indication that the larger input waves (a Q = 0.4 cm, 

Figure 26) peak sooner than the smaller waves. There is also evidence 



*For this case (y < c), photographic measurements indicate that the mode 
is preserved at least over the lower and middle slope areas. 



\ 
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that the measured values fall off to lower values than the theoretical values 
to the left of the zone of peaking (i. e., toward the corner). 

The high-frequency waves over the steeper slope (Figure 28) con- 
form closely to the predicted behavior for x c / x 0 > 0.35, but nearer the 
corner the amplitudes of these waves grow less than predicted. This is 
probably a result of the anomalously large growth of the theoretical solu- 
tions near the corner (recall that for y < c Equation (4-6) is singular for 
x = z = 0). The larger input waves (a Q = 0.4 cm, Figure 28) show a 
tendency to peak earlier than predicted. The measured value at 
x c /x q = 0.35 for a Q = 0. 2 cm is suspiciously low. 

Figures 27 and 30 illustrate the behavior of wave amplitude during 
shoaling of the lower frequency waves for a low slope angle (15 degrees). 
The measured values obviously do not achieve the high predicted values 
well upslope (x c /x q < 0.4 in Figure 27 and x c / x 0 < 0. 5 in Figure 30). 

The departure from theoretical estimates is most noticeable for the very 
low frequencies (co = 0. 047 cps and co = 0. 043 cps) in Figure 27. This sug- 
gests that the waves increase in amplitude as they propagate upslope until 
at some point near the half-distance of the horizontal projection of the slope 
the rate of growth decreases. It was indicated earlier that these waves 
tend to surge as they approach the corner. Possibly a transformation from 
the sinusoidal wave form to a different wave structure (such as a solitary 
wave) that is continually losing energy to viscous dissipation along the slope 
accounts for the nature of the experimental data for these kinds of waves. 

In any event, the amplitude of these lower frequency waves is not as large 
as predicted, and there seems to be evidence of a more gradual dissipa- 
tion of these low frequency waves with the formation of runup. 
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Critical and Supercritical Cases 

The experimental runs for the critical and supercritical cases are 
indicated by a dash in column 6, Table B2 (Appendix B). A comparison of the 
input wave frequency with the critical frequency for each of these runs shows 
that u> c , a necessary condition for this regime. A brief description of 
the reflection properties that define this regime was presented in the first 
section of this chapter. In particular, for case C (supercritical) two fea- 
tures are suggested by the reflection properties: (1) the presence of 

back reflected wave energy; (2) concentration of this back-reflected energy 
in a zone bounded by the bottom slope and the characteristic that intersects 
the corner. Figure 31 further illustrates these features. This figure 
shows that wave characteristics which are incident to the slope from above, 
i. e., those with a downward-pointing vertical component, upon reflection 
from the slope will remain below the characteristic (defined by z = -cx) 
that intersects the corner. The energy flux associated with these charac- 
teristics enters the slope region in a relatively broad band and exits in a 
confined, triangular zone (for example, the triangular area AOB in Figure 
17). As Figure 31 indicates, there is a smaller influx of energy which, 
upon reflection from the flat bottom, propagates upslope within the narrow 
bottom zone (for example, the path of the characteristic "b” in triangle 
AOB). The characteristics associated with this latter energy flux have 
upward-pointing vertical components as they enter their slope regions; 
upon reflection they preserve the sense of the upward-directed vertical 
components. In both cases (characteristics M a" and "b" in Figure 31) once 
the characteristics reflect from the slope the horizontal direction of energy 
is away from the corner. 

Critical conditions (case B) were difficult to achieve experimentally 
for two reasons: (1) y - c required that the wave-maker frequency be 

tuned precisely; (2) it was difficult to determine the Brunt- Vaisala fre- 
quency accurately during an experiment. In addition, the experimental 
groups which involved the 15 degree slope angle had critical frequencies 
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FIG. 31. Reflection of characteristics for case C, y > c where y = tan a, 
z c = |z|, x c = |x| and characteristics are labeled a, b, c. Char- 
acteristic c is "critical” i. e., passes directly to the corner at 
O in (1) and at F in (2). The increased concentration of charac- 
teristics in the triangular zones (hatched areas) is illustrated by 
the path of characteristic a. 
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that corresponded to long horizontal wave lengths over the flat bottom (on 
the order of the tank length). Data for only six runs at the critical fre- 
quency are considered reliable. During each of these runs the input wave 
frequency was held to within two percent of the calculated critical frequency. 

Table B5 (Appendix B) summarizes the experimental results from 
the conductivity-probe measurements for this regime (y > c). Like the 
organization in Table B3, the frequency, period, and probe coordinates 
(x , z )* are given for each experimental run. The experimental runs 

L v 

having the same N and y are again collected into groups; as before, the 
experimental parameters for each experimental group are listed in Table 
Bl. An asterisk in column 2 of Table B5 again indicates that the experi- 
mental run was repeated at least once. The local depth h is shown in 

X 

column 6 of Table B5, and the parameter — „ ~ in column 10 is 

x^ - (z c r 

analogous to the parameter in column 7 of Table B3. 

The experimental method of determining wave number and wave am- 
plitude was discussed in Chapter 3. These quantities are plotted as functions 
of x c /x q and z c /cx c in this section to examine their dependence on (1) hori- 
zontal distance from the slope at a fixed depth, and (2) nearness to z^ = cx c 
(i. e., distance to the theoretical position of the critical characteristic). 

The pertinent theoretical curves were computed from Equations (4-18) and 
(4-20) and are presented in several of the figures for comparison with the 
data. 

It must be pointed out here that qualitative considerations of the 
linear, inviscid solution for y >c, given by Equation (4-16) suggest a priori 
that these solutions might fail to fit the experimental data. It can be shown 
that the velocity components derived from Equation (4-16) become anoma- 
lously large for large coordinate values (x, z). This behavior was not 



In this discussion and in Table B5 x = jx| and z = ) z| for x, z as shown 
in Figure 15. c c 
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observed in the experiments. In addition it is shown later that the theory 
(in particular, Equation (4-17)) selects |n| = 2 as the lowest nonsingular 
inviscid mode permissible over the slope region for each of the experimental 
conditions tested. In the particular experimental geometry, it was not ap- 
parent that the second mode was present over the slope; it was found most 
often that the modal structure was approximately that of the input mode 
( | n | = 1). In spite of these qualitative objections, the following discussion 
of the experimental results is done quantitatively by comparison with the 
theoretical predictions and by consideration of the reflection properties 
(Figure 31). It is felt that this method of presentation aids the interpreta- 
tion of the results. 

Wave number. Figures 32 and 33 illustrate the complicated behavior of the 
experimentally determined wave numbers for this case. Data from three 
experimental groups and for frequencies co = 0.051 cps and co= 0.066 cps 
are represented in these figures. The calculated positions of the critical 
characteristic (determined by z^ = cx ) for two of the experimental groups 
are shown along the horizontal axes. The two theoretical curves shown in 
each figure were computed for two different values of the mode number (see 
qualitative discussion at end of previous section). Equation (4-20) was used 
to determine the minimum permissible nonsingular mode number for the 
curves labeled n = 2. The mode number n = 1 was preselected for the theo- 
retical curves labeled n = 1. Although this latter choice violated the con- 
dition necessary for finite values of the velocity components and their first 
derivatives at z = -cx (Equation (4-20)), it nevertheless forces the theoretical 
mode number over the slope to be equal to the input mode of the experiments. 

It is apparent in the figures that there is no straightforward relation- 
ship between linear theory (n = 2) and the data. There is only qualitative 
indication that the curves for n = 1 have the same trend as the experimental 
results for x c / x 0 — 0. 035 in Figure 32 and for x c /x q < 0. 45 in Figure 33. 

It is obvious that the theoretical estimates for n = 1 grow toward very large 
values at the critical points (for example, X C / X Q = 0. 26, Figure 33). The 
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Fig. 32. 

Upslope variation of wave number, y > c. 
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Fig. 33. 

Upslope variation of wave number, y > c. 
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experimental values near the calculated critical points are apparently 
bounded. The dip in the data in Figure 32 for X J X Q - 0-30 is not apparent 
in Figure 33. Possibly the absence of the dip in the data from experimental 
group 4 (Figure 33) can be explained by considering nonlinear effects. The 
measured value of e m (Table B5) for this particular data set is approximately 
0. 5 at x c / x 0 = 0.36. This suggests that nonlinear amplitude growth might 
be significant near these positions. By comparison, the measured value of 
e m for the data in Figure 32 is only 0. 13 for X C / X Q ~0. 2 (Table B5). It 
should also be noted that there is a significant shift in the peak value to the 
right of the calculated critical point in Figure 32. Possibly this shift indi- 
cates that the singularity is removed by nonlinear or viscous processes and 
that the linear inviscid solutions are not valid in this particular physical 
situation. 

Data from experimental groups 3 and 4 are plotted on log-log scales 
in Figures 34 and 35, respectively. Theoretical curves for n = 1 are again 
shown in these figures. The qualitatively similar trend between the theo- 
retical curve and the data is apparent, particularly in Figure 35. Values of 
wave number for the larger amplitude input waves (a Q = 0.4 cm) tend to be 
higher than those for the smaller waves (a Q = 0.2 cm) in experimental 
group 3 (Figure 34). 

Figure 36 shows the results of experimental group 3 plotted as a 

function of z /cx . It is obvious that the theoretical curve for n = 1 grows 
c ' c 

toward infinite values of k /k at z = cx : the theoretical curve for n = 2 

m o c c 

(obtained from Equation (4-17)) is finite at this point. The trend in the 
data values reverses at about z c /cx c =* 0.80, and again at z c /cx c ^ 1.0. 

No simple relationship between the data and the theoretical curves is 
apparent from this figure. 

Figure 37 summarizes the experimental results of wave number for 
case B (y = c). The log-log plot emphasizes the steep rise toward larger 
values of wave number for decreasing distance to the slope. It is interesting 
that the magnitudes of wave number for the larger input waves (a Q = 0.4 cm) 
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Fig. 34. 

Upslope variation of wave number on log-log plot, y > c. 
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Fig. 35. 

Upslope variation of wave number on log-log plot. 
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Fig . 36 . 

Have number versus magnitude of z c /cx c , y > c 
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in experimental group 3 are larger than the corresponding values of the 

smaller input waves for x /x <0.3, but for x /x >0.3 the wave numbers 

c o c o 

of the larger waves are consistently less than those of the smaller waves. 

In general, the wave number variations with distance to the slope for this 
case are similar to the variations shown in Figure 34 for the supercritical 
case. No theoretical comparison is shown, since for y = c the minimum 
permissible nonsingular inviscid mode over the slope, given by Equation 
(4 -17), is infinite. 

In summary, the results show a complicated relationship between 
wave number and position over the slope. There is evidence in Figures 
34 and 35 that the wave number increases rapidly as the waves pass over 
the lower and middle sections of the slope for the supercritical case, but 
there does not appear to be a systematic pattern to the data near z = cx . 

The variations in wave number for the critical case are qualitatively simi- 
lar to those for the supercritical case; a large increase in the wave number 

for smaller x /x is apparent. The trends of the theoretical values of wave 
c' o 

number computed for n = 1 are similar to the measured values for z c /cx c < 1. 

It is obvious that the unbounded theoretical values at z = cx for n = 1 are 

c c 

removed when the condition given by Equation (4-17) is satisfied. For both 
cases, the measured values of e m (Table B5) are greater than 0. 1 for 
most of the measurements taken closest to the corner. This suggests that 
nonlinear effects might be important near these positions. 

Wave amplitude. The experimental estimates of wave amplitude were ob- 
tained from the periodograms of the conductivity-probe measurements 
(Chapter 3). The results are presented as ratios of the square of the wave 
amplitude over the slope to the square of the wave amplitude over the flat 
bottom for various frequencies of the input waves. Like the earlier presen- 
tation for wave number, these results are illustrated graphically with 

z /cx and x /x used as the horizontal axes, 
c c c' o 
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The results for wave amplitude are easier to interpret in terms of 

the reflection properties of the slope than were the results for wave number. 

Earlier theoretical discussion indicated that near the characteristic defined 

by z = cx and in a zone between this characteristic and the bottom slope 
c c 

(Figure 31) there is a possible region of concentrated energy flux. Figure 

38 shows that the measured values of wave amplitude are highest near 

z /cx =1.0 for each frequency tested in this group except w= 0.039 cps. 

However, it is not conclusively shown in this figure that the amplitudes peak 

exactly at z = cx . This figure also shows that the ordinate values for the 
c c 

various frequencies at equal values of z /cx generally increase with de- 

creasing frequency. In the light of Figure 31, back-reflected wave energy m 

might contribute smaller amounts to the total energy at locations above the 

critical characteristic (i. e., at positions given by z < cx ) for frequencies 

c c 

nearer to critical. This in turn might account for the lower measured 
2 

levels of (a /a Q ) f° r the frequencies nearer to critical in Figure 38. 

Figure 39 is a similar plot for data from experimental groups 3 and 4; 

two input wave amplitudes are represented in the values from experimental 

2 

group 3. Also shown in Figure 39 are theoretical curves for (p /p ) de- 

w o 

rived from Equations (4-10) and (4-18) by the methods described earlier. 

The mode number |n| was determined from the condition in Equation (4-17) 

for curve n = 2 in this figure. The theoretical values (a /a ) decrease 

m o 

toward the slope and undergo a smooth transition through the critical point 
z /cx = 1.0, where they then begin to increase toward the slope. The 
theoretical curve with n = 1 has increasing values toward the slope until the 
critical point (z = cx ) is approached; a sharp dip in the theoretical values 
at this point followed by a rapid rise suggests that the solutions are not well 
behaved at this critical point. The ordinate values for the two input ampli- 
tudes of group 3 in Figure 39 agree within 10 percent of one another for 
values of z /cx < 1.0; in general, the larger input waves have slightly 
higher measured amplitudes in this range. It is also apparent that the 

larger input waves grow more rapidly as the value z /cx = 1.0 is approached. 

c c 

Owing to the larger values of e m (Table B5) for the waves of larger input 
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Fig. 38. 

Wave amplitude versus magnitude of z c /cx c , y > c. 



91 




cd = 0.051 c ps 


EXP group 


0 o , cm 


3 


O 0.2 
A 0.4 


4 


□ 0.2 



Fig. 39. 

Wave amplitude versus magnitude of z c /cx c , y > c. 



92 



amplitude (a Q = 0.4 cm), the increased disparity in the measured amplitude 

values for input waves of 0.4 and 0. 2 cm, respectively, at values of 

z /cx >1.0 in Figure 39 might be the result of nonlinear effects. 

0 0 

Figure 40 shows the experimental results and theoretical curves for 

co = 0. 051 cps (group 3) on a log-log scale; the horizontal axis is distance x 

to the corner normalized by the total horizontal projection of the slope x q . 

The rapid increase in amplitude upslope for 0. 28 < x /x < 0. 60 is readily 

c o 

apparent in the data; however, the theoretical values for n = 2 show an op- 
posite trend, toward lower values for decreasing distance to the slope. The 
measured amplitudes below the indicated position of the critical point remain 
relatively high; this suggests amplification of the motion near the bottom. 
Nonlinear effects also might be important near the corner (x c /x q = 0. 19) 

The results for the probe measurements of wave amplitude at critical 

frequencies are shown in Figure 41. The increase in amplitude over the 

range 0. 28 < x /x < 0. 60 is similar to the increase shown in Figure 38 for 
c o 

the same range. However, the maximum amplitude measured for these 
waves is considerably less than the corresponding value for the supercritical 
case (Figure 40). Larger viscous dissipation along the slope in the critical 
case might explain the smaller peak amplitude levels. 

In summary, the data indicate that the wave amplitudes generally 

increase toward the slope; the higher frequency waves (i. e., those closer 

to u> c ) for these cases show rapid growth during shoaling over the middle 

sections of the slope. Peak levels of amplitude were generally found near 

z /cx ta 1.0, with lower measured peak values for waves of critical fre- 
c c 

quency. In both cases the motion is relatively large near the slope. In 
general, the conductivity-probe measurements of amplitude agree qualita- 
tively with the reflection properties of internal waves from a rigid, sloping 

bottom. There is some indication from the measured values of e that 

m 

nonlinear amplification of the amplitude occurs for small values of x c /x q . 
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Upslope variation of wave amplitude, y > c. 
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Upslope variation of wave amplitude, y > c. 
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It might be added here that the bounded experimental values at 
z = cx suggest that viscous dissipation and nonlinear processes might act 
to remove the singularity at the critical characteristic. For example, the 
shift in the measured peak wave number in Figure 32, mentioned earlier in 
this section, and the formation of vortices near the bottom for these cases, 
discussed in the next section, suggest that the linear inviscid model is not 
properly tested by these experimental results. Nonlinear and viscous effects 
in the interior, particularly in the vicinity of the critical characteristic, 
should be examined in future analytical treatments of this problem. 

BOUNDARY ACTIVITY 

The discussion of the bottom boundary layer and the zone of breaking 
follows the procedure of separating the experimental results into two hydro- 
dynamic regimes (Iy < c and II y > c). The distinction between these two 
regimes was defined in Figure 14. Qualitative features observed during the 
experiments and illustrated by several photographs are presented at the 
outset as an introduction to the later discussion of quantitative results. In 
the discussion of quantitative results, experimentally determined velocities 
in the bottom boundary layer at several positions along the slope are shown 
to agree with theoretical values computed from the linear solutions in 
Equation (4-24) for case A. On the basis of the good agreement between 
experiment and theory for the conditions tested in case A, Equation (4-24) 
is used in Chapter 5 to develop a model for sediment movement induced by 
shoaling internal waves. 

Case A: y <c 

Qualitative results. The shadowgraph images and the motion of neutrally 
buoyant floats permitted excellent real-time observations of the oscillatory 
flow in the boundary layer and breaking zone. One obvious major feature 
of this flow was its intensification near the bottom along the slope as com- 
pared with the flow along the horizontal bottom. The maximum velocities 
generally increased at positions farther upslope, with a sharp growth in the 
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velocities near the zone of breaking. The oscillatory flow remained laminar 
in the boundary layer up to the zone of breaking. 

The breaking of waves on the relatively steep slopes (30 and 45 

degrees) is characterized by considerable turbulence and vertical motion 

below the free surface and by the generation of spatially irregular fine 

structure in the density field. This kind of breaking is also observed for 

relatively high-frequency waves on the smaller slopes (7 and 15 degrees). 

In contrast, the lower frequency waves (i. e., those waves with frequencies 

closer to co ) break much less violently, and typically form surges that run 
c 

up the slope for relatively long distances. 

More detailed observations of the higher frequency waves indicate 
that shortly after the first arrival of the waves near the top of the smaller 
slopes (a = 7 and 15 degrees), a small vortex forms along the slope be- 
neath the crest of the advancing wave form. Figures 42 and 43 illustrate 
the onset of this instability beneath the wave crest. In Figure 43 the maxi- 
mum amplitude of the incoming wave at the entrance to the slope region 
was approximately 0.4 cm, whereas the corresponding amplitude in Fig- 
ure 42 was 0. 2 cm. The time that had elapsed from the start of the wave 
maker is shown in the guide sheet preceding Figure 42. Subsequent devel- 
opment of the instability in Figure 42 is shown in Figures 44a and 44b, and 
later stages of the instability in Figure 43 are shown in Figure 45a and 45b. 
The lamellar structures in these photographs represent small-scale depar- 
tures from the constant vertical index-of-refraction gradient characteristic 
of the stratified fluid with no wave motion. A constant index-of-refraction 
gradient (or equivalently, constant density gradient) produces a uniform 
uninteresting shadowgraph image (see center of field in Figure 42). In 
spite of the apparently chaotic density structure, the basic wave-form can 
still be discerned in the photographs. Qualitatively, the larger amplitude 
waves show a higher intensity of breaking, as evidenced by the contrast be- 
tween Figure 44b and Figure 45b. Not shown in this series of photographs 
is the protrusion of thin horizontal layers back into the interior region 
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over the slope. The protrusions lose their intensity (as measured by 
photographic contrast) over very short horizontal distances (5 to 10 cm) 
and are confined to the upper vertical levels, where the breaking occurs. 

The protrusions, thin tongues of mixed fluid that are weakly converted back 
into the interior, are probably caused by mixing in the zone of breaking. 

In contrast to the higher frequency waves, the waves of lower fre- 
quency on a small slope (7 degrees or 15 degrees) show less mixing and 
fewer laminae in the zone of breaking. Each wave forms a surge in which 
the velocities of neutrally buoyant particles within the surge become approxi- 
mately equal to the wave celerity as the wave progresses upslope, somewhat 
analogous to solitary surface waves. The surging motion eventually produces 
a runup and backwash zone along the bottom very near the corner. The runup 
appears as a thin streamer in the shadowgraph image and represents a pro- 
trusion of fluid of higher density into the near- surface layer. Figures 46a 
through 46c illustrate the development of this phenomenon for waves whose 
initial amplitude is 0.4 cm, Figures 47a through 47c show waves whose 
initial amplitude is 0. 2 cm. Vortex instabilities initially form beneath the 
wave crests and propagate upslope with the wave form. The quasi-steady 
development for the 0. 4 cm and 0. 2 cm waves is shown in Figures 46c and 
47c, respectively. In these figures, the surging motion essentially fills 
the water column in the left portion of the photographs. Farther upslope a 
vortex forms and eventually dissipates into runup. 

Figures 48a through 49c show a corresponding development of break- 
ing for waves of 10. 7 sec period and 0.4 cm initial amplitude over a 30- 
degree slope. The mixing becomes very intense after the collapse of the 
initial vortex (Figure 48a). It is interesting to note the reversal in the cir- 
culation of the lead vortex as it propagates upslope and collapses. The re- 
versal is a consequence of the shear between the surging motion near the 
bottom and the opposing flow above the vortex. The following sketch in 
Figure 50 illustrates this pattern. 



Figs. 42 through 49. Guide to photographs of shadowgraph images in the 

following figures. 
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Fig. 42. 




Fig. 43. 
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Fig. 44 a. 




Fig. 44 b. 
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Fig. 45 a. 




Fig. 45 b. 
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Fig. 46 a. 




Fig. 46 b. 




Fig. 46 c. 
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Fig. 47 a. 




Fig. 47 b. 




Fig. 47 c. 
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Fig. 48 a. 




Fig. 48 b. 
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Fig. 49 a. 




Fig. 49 c. 
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FIG. 50. Advancing initial vortex at successive locations A, B, and C. 
Note observed reversal in circulation at position C. 




The outflow above the vortex balances the flow of water that is advected 
toward the corner along the bottom by the wave motion. 

As indicated by the well developed lamellar structures in Figures 
49b and 49c, the steeper slopes appear to produce a high intensity of breaking. 
The formation of a surge with associated runup was not observed for the 
steeper slopes (30 degrees and 45 degrees) for any wave frequency tested. 

Figure 51 presents a summary of the observed positions of breaking 
on slopes of 30 degrees. The accuracy of the distance is low, owing to 
the difficulty in defining the point of breaking; however, the relative order 
of the values of for the waves of different frequencies shown in the figure 
is reliable. The general trend indicates that for a given N, the higher fre- 
quency waves break first, i. e., in deeper water. 
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FIG. 51. 
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Quantitative results. The experimental techniques that were used to mea- 
sure the motion in the boundary layer are discussed in Chapter 2. The 
results of these measurements are compared with the predictions of the 
linear, viscous boundary layer solutions in order to check the validity of 
these solutions for certain conditions and provide some empirical estimates 
of the motion where the theory does not apply. 

Since successive displacements of the dye streaks and neutrally buoy- 
ant particles were measured at one-second time intervals, the velocity 
values computed from these measurements are essentially one- second time 
averages in the Lagrangian sense. Equation (4-24) provides a linear, 
Eulerian solution for the velocity field parallel to the bottom in the viscous 
boundary layer; this solution is identical with the Lagrangian one to first 
order in e. The Lagrangian velocity u^ can be determined in principle 
from the Eulerian velocity u^,: 

u L (a, t) = u E (a,t) + J^ Q u(a, t') dt' • V- u E (4-27) 

where a is a position of the dye streak, say, at timetb; t - t Q is small com- 
pared to the wave period T, and V is the gradient operator (Phillips, 1966). 
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This expression for u T is correct to 0(e) if u_ is correct to this order. 

Li Li 



Let 



- - (1) 2 ~* (2) 
U E = €U E + 6 U E 



(4-28) 



where u^, is the linear solution for the velocity field. It is easy to show 
that there are two second order contributions to u^ in Equation (4-27): 

(1) an d (2) ( Ugdt') * Vu^,. The first contribution (1) is due to the 

nonlinear terms in the equations of motion; the second (2) is commonly re- 
ferred to as the Stokes velocity and can be computed in principal from Equa- 

o 

tion (4-21). Hogg and Wunsch (1970) have shown that to order e the time 

average of (T for this problem is zero. This implies that there is no net 
L 

transport of water particles induced by internal waves near a sloping boun- 
dary. It also implies that since u^(^) = 0, then 

u = -( u E ^W) • Vu E ^ (4-29) 

o 

where the bar denotes a time average over one wave period. 

A complete analysis of the data was carried out for three wave fre- 
quencies in regime I (y<c). Two of these analyses involved waves on a 

30 degree slope (o> = 0. 105 cps and w = 0.094 cps, group 7, Figure Bl, 
Appendix B) and one involved waves on a 15 degree slope (co = 0. 105 cps, 
group 2). For the sake of brevity, and without any appreciable loss of de- 
tail, only the results from co = 0. 105 cps will be discussed at length. 
Instances where measurements from the other runs showed significant 
features of deviations from this one will be mentioned in the discussion. 

The results for regime H ( y > c ) are discussed in a later section. 
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A time- sequenced series of streak and particle positions from 
which velocity information was obtained is shown in Figure 12. The streaks 
are numbered according to consecutive, one-second time marks relative to 
the initial one. The numbers on the particles correspond to those for the 
streaks. The increase in displacement of the streaks nearer to the bottom 
agrees qualitatively with the amplification predicted by Wunsch (1969) for 
y <c (Figure 52). 

Figure 53 summarizes some of the relevant parameters that were 
estimated from the measurements. The five station marks in the diagram 
show the locations of dye pellets. The coordinates for these locations are 
given in columns 2 and 3. The definitions of the various Reynolds numbers 
and estimates of boundary-layer thickness are given in the legend. The 
approximate maximum error in the measured estimates of the boundary- 
layer thickness (<5 ) is 10 percent of the listed value. 

The boundary-layer behavior for steep slopes (30 degrees) is anom- 
alous on the lower portion of the slope. The data at station 4 for = 0. 105 
cps (Figure 57) and go = 0.094 cps both reveal erratic streak displacement. 
The estimates of boundary-layer thickness at this location were slightly 
larger or approximately equal to those at station 3. The measured thick- 
ness of the boundary layer at corresponding positions on the smaller slope 
(15 degrees) generally increased in an upslope direction to the zone of 
breaking. The anomalous behavior for steeper slopes was probably associ- 
ated with an adjustment in the modal structure as the waves progressed 
from the flat-bottom region to one of appreciable slope. The irregularities 
in the motion at station 4 are discussed later in this section. 



The larger values of 6 as compared with 6 ^ at stations 1, 2, and 

3 represent significant departures of the measured estimates from the 

theoretical values. In addition, the theoretical value is constant over the 

entire slope, whereas the measurements indicate that <5 m varies slightly 

with slope positions. Over the flat bottom, 6 is slightly larger than the 

1/2 m 

value given by (2y/co) ' , although this difference was not significant for 
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Fig. 52. Theoretical prediction of increase in 
velocity components along the slope, 
are arbitrary (after Wunsch, 1969) . 
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FIG. 53. Estimated of boundary layer thickness and 
Reynolds numbers. r < c. 
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the experiments in which co = 0. 105 cps, a - 15 degrees. It appears that 
the boundary layer thickness given by linear theory gives the correct 
order, mm rather than cm, but the measured values tend to be significantly 
higher. However, it must be recalled that the definition of boundary layer 
thickness for the purposes of measurement, i. e., the distance between the 
bottom and the point of inflection in the velocity profiles, is somewhat 
arbitrary, and might account for the discrepancies with the theoretical 
values. 



It should also be mentioned here that at station 1, the computed 
value of e is relatively high (g ^ 0. 5); nonlinear effects are probably im- 
portant there. It was mentioned earlier that Wunsch (1970) recently showed 
theoretically that there can be a mean upwelling velocity independent of wave 

motion in a layer whose thickness is 5 : 

J w 



w 



<„ k) 1/4 /n 1/2 



where v and k are the vertical diffusivities of momentum and heat, re- 
spectively. In terms of molecular processes, v ~ lo2 cm ^/ sec and 
k ~ 10 - 3 cmV sec . Using the value of N given in Figure 53, N ~ 1 sec \ 
we can estimate 8 ~ 0.06 cm for this particular experiment. 

The curves in Figures 54 through 58 show velocity versus depth for 
all five stations. Each curve represents a one-second time average of the 
Lagrangian velocity parallel to the bottom. The consecutive numbers on the 
curves indicate the time history of this velocity over one wave cycle. The 
curves were computed only to the vertical level nearest the bottom for 
which the streak data was considered reliable. The horizontal velocity 
scales are uniform in each of the figures except for the expanded scale of 
the flat bottom results (Figure 58) and the condensed scale at station 1; the 
vertical scales are the same for all figures. 

It is important to note that at each of these stations the observed 
motion in and near the boundary layer was essentially parallel to the bottom. 



UJ o 
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Fig. 54. Vertical distribution of velocity near the bottom at station 
a Q = 0.2 cm, OJ = 0.105 cps , N = 0.155 cps , oL =30 deg, 



114 



in 




Fig. 55. Vertical distribution of velocity near the bottom at station 
a 0 = 0.2 cm, CO = 0.105 cps , N = 0.155 cps , OC. =30 deg, 
x = 23 . 0 cm. 
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Fig. 56. Vertical distribution of velocity near the bottom at station 3 
a Q = 0.2 cm, QJ = 0.105 cps , N = 0.155 cps , OC =30 deg, 
x = 36.0 cm . 
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Fig. 57. Vertical distribution of velocity near the bottom at station 4 
a 0 = 0.2 cm, CO = 0.105 cps , N = 0.155 cps , c<, =30 deg, 

x = 51.0 cm. 
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Only very small velocities normal to the bottom were detected from the 
motion of the plastic beads, and this vertical motion was usually confined 
to regions adjacent to and in the zone of breaking. 



A reason for the anomalous behavior of the velocities at station 4 
(Figure 57) was hypothesized earlier in this section. An alternative ex- 
planation is that the abrupt change in bottom geometry acts as a source of 
diffracted waves which locally affect the pattern of motion. However, 
since the density did not vary linearly with depth at the bottom (for example, 
Figure 9), but instead gradually became more homogeneous in a small 
layer near the bottom, the dynamics in this transition region are probably 
quite complicated. Visual observations of the motion of neutrally buoyant 
particles and dye patches that were placed at the transition in bottom geom- 
etry during several experiments (for 15 and 30 degree slope angles) did not 
disclose any significant patterns other than those already mentioned. 

Figure 59 shows a comparison between the theoretical and experi- 
mental RMS values of the velocities at selected vertical levels shown in 
Figures 54 through 58. The experimental RMS value at each level was 
computed as: 



U RMS 





where 



u 

MS 






i = 1, 2 T 



and T was chosen as the smallest integer greater than the wave period. 

For this particular set of data, T = 10 (wave period ^ 9.5 sec). For 
stations on the slope, theoretical RMS values were obtained from the Eulerian 
expression for the boundary layer velocity component parallel to the slope 
(Equation (4-24)): 

u RMS (the0ry) = (1 /2 “ e ' “ E ) 1/2 ' 
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These RMS values are listed in Figure 59 as first order theoretical estimates 
of the RMS Lagrangian velocities. 

Qualitatively, we might expect that the measured values are more 
accurately represented by the sum of the first order Eulerian values and 
the second order values consisting of the Stokes velocity component and the 
nonlinear component, particularly farther upslope where e become larger. 

When the various sources of uncertainty in the measurements of 
velocity and experimental parameters (A q , N, co, h, etc.) are considered, 
a rather good comparison would require that the theoretical and experimental 
values fall within 15 percent of one another at each level ri . The results in 
Figure 59 show that excluding station 4 there is excellent agreement (with- 
in 10 percent) between the measured RMS values and the corresponding 
first order, theoretical values. The excellent comparison at station 1 is 
particularly remarkable, considering that the estimated value of e is approx- 
imately 0. 5 there. The high measured values at station 4 are a consequence 
of the anomalous motion that has already been discussed. 

The results for the maximum velocity |u| max parallel to the bottom 
at the selected levels (columns 6, 7, 8, Figure 59) also show good agreement 
between experiment and theory. Theoretical estimates of the maximum 
values of this velocity component were again obtained from Equation (4-24). 

The experimental values are generally higher than the theoretical first order 
values. The large deviations between experiment and theory at station 4 are 
again obvious. The results for the other analyses of this kind show similar good 



FIG. 59 
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agreement between theory and experiment. The velocity values for station 
4 with the 15 degree slope are more in agreement with the theoretical 
values (within 25 percent) than are the results in Figure 59 and did not show 
anomalously large velocity values for this station. 

The asymmetry in the velocity profiles at station 1 (Figure 54) 
contrasts with the more symmetrical curves for the other stations. 

Station 1 is located well upslope and only a few centimeters from the zone 
of breaking. The skewness toward positive velocities suggests that there 
might be a measurable net upslope motion of the water particles over one 
wave cycle. Although the maximum velocity is directed downslope (curve 
1, Figure 54), the magnitudes of the upslope component are quite high 
(approximately 1 cm/ sec) for a considerable fraction of the wave cycle. 

The local measured value of € is about 0. 5 at this station, which suggests 
that nonlinear effects might be important. It is also interesting to note 
that S m at this station shows a decrease from that at station 2 (Figure 53). 

Experimental estimates of net drift at the outer edge of the boundary 
layer were computed at each station by algebraically averaging the mea- 
sured velocity values at the estimated level of 6 m in Figures 54 through 58, 
respectively. The average is for one wave cycle (sampling interval of one 
second) and represents a single estimate of mean streaming at each station. 
No statistical significance can be inferred from this method, but the mea- 
sured results in Figure 60 show that the streaming is insignificant at all 
slope stations except station 1. Here the net drift is approximately 1. 3 mm/ 
sec. This drift is probably induced by mixing in the nearby (3-5 cm) zone 
of breaking. As discussed earlier, protrusions of thin fluid laminae slowly 
advect into the interior over the slope from the zone of breaking. The direc- 
tion of the drift measured at station 1 is toward the breaker, possibly indi- 
cating that this flow compensates in part for the mass efflux in the protru- 
sions. The observed direction net drift near the free surface above this 
station is also toward the zone of breaking. 
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FIG. 60. Measured values of mean streaming at 15 = 6 
stations 1 through 4 on slope. 



Station No. 


x , cm 
c’ 


I, cm 


u g , cm/s< 


1 


8,7 


10.05 


+0. 13 


2 


23.0 


26.6 


-0.07 


3 


36.2 


41.8 


-0.003 


4 


51.0 


58.9 


-0.023 
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Owing to the poor reliability of the raw streak data very near the 

wall (<0.07 cm), the best experimental estimate of bed shear stress, T , 

’ o’ 

was obtained by extending the reliable portion of the velocity-depth curves 
to the zero point at the bottom. The shear stress was then computed from 
a simple straight-line fit between the zero point and the lowest point on the 
velocity curve. This gives a low estimate of the bed shear stress, since 
the actual shape of the velocity curves in the boundary layer is more 
exponential. The maximum bed shear stress /T Q /max measured at each 
location is listed in column 7 of Figure 59; the time variations in t q over 
one wave cycle corresponding to the time variations in the velocities in 
Figures 54 through 58 are shown in Figure 61(a-d). The theoretical values 
of bed shear stress were computed from Equation (4-24): 

T 0 - M du/drf , rf = 0 

where 

- = 0.01 — — — . 

cm-sec 

The maximum theoretical value is shown in column 8 , Figure 59, 
for comparison with the corresponding measured value. Since the theo- 
retical estimates were made from the linear Eulerian expression of velocity 
at each station, they correspond to the experimental Lagrangian estimates 
only to first order in e. It can be shown that n varies by only 4 percent over 
the density range for the particular experimental conditions shown in Fig- 
ure 53. In spite of the poor comparison of the estimates of boundary layer 
thickness given by linear theory and experiment, the theoretical estimates 
of maximum bed shear stress at stations 1-3 compare reasonably well with 
the measured values. This is apparently the consequence of the higher ex- 
perimental values of maximum velocity as compared to the first-order 
theoretical values for levels near the boundary. Figure 61 also shows the 
excellent fit of the theoretical curve computed from first order solutions to 
the experimental values at stations 1, 2, and 3. 



124 



Shear stress, dynes/cm 2 




. 61 a. Experimental and theoretical values of bed 
shear stress x 0 at station 1. The experi- 
mental values were obtained for each 
second of time t during one wave cycle. 



Fig 
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Shear stress, dynes/ cm 2 




Fig. 61 b. Experimental and theoretical values of bed 
shear stress t q at station 2. The experi- 
mental values were obtained for each second 
of time t during one wave cycle. 
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Shear stress, dynes /cm 2 




Fig. 61 c. Experimental and theoretical values of bed 
shear stress t q at station 3. The experi- 
mental values were obtained for each second 
of time t during one wave cycle. 
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2 

Shear stress , dynes /cm 




Fig. 61 d. Experimental and theoretical values of bed 
shear stress t q at station 4. The experi- 
mental values were obtained for each second 
of time t during one wave cycle. 
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Cases B and C: y = c and y > c 

The description of the boundary- layer motion in this section follows 
a format similar to that presented in the discussion of boundary activity for 
case A (y < c). The qualitative remarks deal mainly with the observations 
and photographs of shadowgraph images, whereas the quantitative results 
were derived from the measurements of the motion of dye streaks and 
neutrally buoyant particles. 

Qualitative During those runs for the critical case (y = c) with moder- 
ately steep slopes (30 and 45 degrees), the zone of breaking near the cor- 
ner was characterized by intensive mixing as evident in the shadowgraph 
images, similar to that described for the high-frequency waves on a 30 
degree slope (see earlier discussion for y <c in this chapter). For runs 
with a slope of 15 degrees, each wave, for both the critical and supercriti- 
cal cases, formed a surge that usually ended with a short runup. This 
runup usually occurred at the bottom of the more homogeneous near-surface 
layer, if present, and appeared like a gentle, lapping motion with much 
smaller upslope penetration than the runup for the subcritical case, which 
persisted over a much longer distance. In summary, the intensity of 
breaking and runup for y > c is decreased as compared with the experi- 
mental observations of these processes for the subcritical case. 

Figures 62 through 64 illustrate in a very striking way that the in- 
ternal wave motion for these cases (y 2: c) is associated with the growth of 
an instability along the slope. The shadowgraph observations of the insta- 
bility illustrated in these figures indicate that a line of regularly spaced 
vortices forms near the bottom. These vortices grow and decay over each 
half cycle as they oscillate along the slope at the frequency of the input 
waves. It appeared that on the downslope propagation during the wave 
cycle, the vortices would "flatten", indicating that they form only during 
the upslope movement of the water particles. This is indicated in Figures 
63a and 63b. The circulation of adjacent vortices has the same sense; 



Figs. 62 through 64. Guide to photographs of shadowgraph images in the 

following figures. 
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Fig. 62 a. 




Fig. 62 b. 
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Fig. 63 a. 




Fig. 63 b. 
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Fig. 64 a. 




Fig. 64 b. 
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this circulation is indicated diagrammatically below. The axis of each 
vortex is oriented approximately normal to the side walls of the tank, 
and the maximum distance observed between the core of adjacent vortices was 
was about 2 cm. Unlike the vortices which form beneath the shoaling, 
high-frequency waves well upslope and which precede the violent breaking 
(Figure 45), these vortices retain their cyclic growth and decay pattern 
for the duration of the input wave motion. Their regular spacing was much 
shorter than the wave lengths of the input waves, particularly at the middle 
and lower slope regions. 




Shortly after the formation of the line of vortices, thin streamers or 
wisps appear near the cusp of each vortex and penetrate horizontally back 
into the interior of the fluid (note the thin lighter bands that alternate with 
the darker bands in Figures 62 through 64). Qualitatively, it appears that 
fluid of slightly different densities from vertical levels adjacent to each 
vortex cell is mixed in the cell and returns to the interior as a thin streamer 
at a new equilibrium level. The result of this process is alternate layers of 
thin, nearly homogenous streamers and thicker, linearly stratified laminae. 
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Quantitative (1.) Vortices and Layering. The onset of vortex activity for 
any particular run normally occurs within 10 to 20 wave periods following 
the arrival of the first wave motion at the corner region of the slope. The 
initial vortex generally appears somewhere along the middle section of the 
slope, and the instability spreads gradually toward the corners. The ini- 
tiation of vortex activity was observed as a function of input frequency for 
experimental groups 4 and 5; the results are shown in Figure 65. These 
results indicate that, for this set of experiments, (1) the vortices occur 
earliest (in time measured from the start of the wave maker) for co =* co 
even though the lower frequency waves have a higher group velocity; (2) 
there is some co < co c for which the vortices do not appear, and (3) there 
is a small range of co slightly greater than co (i. e., waves in the sub- 
critical range) for which vortices form. 

The diameters of the vortices were difficult to measure accurately 
from the shadowgraph images because of the constant process of either 
growth or decay during each half cycle. From data taken from various 
photographs, the size of the vortex diameters is largest for co =* co Q and 
generally decreases for lower and higher frequencies. Figure 66 shows 
the size variation of maximum vortex diameters measured at midslope 
positions during experimental group 4. The corresponding input wave 
frequencies are listed in column 1. 

The position of the vortices relative to the bottom also varies over 
each half cycle of the wave motion. As a vortex grows in diameter it also 
appears to move away from the boundary; conversely, as its diameter 
diminishes its distance from the bottom decreases. The maximum distance 
measured between the center of a vortex and the bottom was approximately 
1. 2 cm. * 



It is worth noting here that the alternate motion toward and away from the 
boundary by a diminishing and enlarging vortex cell gives the impression 
of a periodic pumping in a direction normal to the bottom. 
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FIG. 65. Onset of vortex motion as a function of input 
wave frequency. 
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-1 

co, sec 


T, sec 


d, cm 


0.546 


11. 5 


0.32 


0.547 


11.7 


0.39 


0.519 


12. 1 


0.43 


0.495 


12.7 


0.52 


0.487 


12.9 


0.51 


0.480 


13. 1 


0.52 


0.468 


13.4 


0.47 


0.458 


13.7 


0.47 


0.430 


14.6 


0.43 



FIG. 66. Maximum size of vortex diameters 
as a function of frequency 



go = 0.495 sec 
c 

N =0.99 sec 
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Shadowgraph images of the thin, horizontal streamers (for example, 
Figure 62b) usually intensify with time if the wave motion continues. Their 
initial appearance generally varies with frequency much in the same way as 
the onset of the vortices. The streamers are usually sharply outlined as 
shadowgraph images near the vortices, but within short horizontal distances 
from the slope (order of 2 - 4 cm) their images weaken. Although they can 
be detected out to distances of 20 cm from the bottom from the shading 
variations on the shadowgraph, the contrast between the images of the 
streamers and adjacent layers becomes very indistinct at distances greater 
than 20 cm from the slope. Their deepest observed penetration of the 
interior was approximately 25 cm, measured horizontally away from the 
slope. This particular measurement was made after one hour of continu- 
ous wave motion at the critical frequency for a run during experimental 
group 4 {a - 30 degrees). The thickness of the streamers during this 
particular run was 2 mm near the boundary, and the vertical distance 
separating two streamers was typically 5 mm. Normally the vertical 
separation between streamers is approximately equal to the diameter of 
the nearest vortex. The shadowgraph images of the streamers are most 
distinct for wave frequencies near a> c > 

Several conductivity probe traverses through short vertical sections 
(3 and 4 cm) were made in order to detect changes in the linear density- 
depth profiles that might be associated with the streamers. Figure 67 shows 
the measured structure for a typical vertical traverse of a probe after the 
formation of the streamers. Small scale irregularities in the curve (after 
the wave motion) are typical of these traverses. The values along the 
horizontal axes represent changes in output voltage relative to the null po- 
sition that was chosen approximately at the center of the traverse; the 
vertical axis is actual tank depth in cm. Changes in equivalent density de- 
termined from the calibration curves are also shown along the horizontal 
axis. Although sharp changes in density at the edges of the streamer are 
apparent from the photographs, some conductivity probe and specific gravity 
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Fig. 67. 

Density-depth profile through a region of 
streamers (numbered) and intervening layers. 
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measurements suggest that the density is not homogeneous within the 
streamer, but has a small, nonzero gradient there. This latter remark is 
based on only six measurements within two streamers (not shown in Figure 
67), and these were taken near the associated vortex cell close to the slope, 
where the streamer could easily be identified from shadowgraph observa- 
tions. The apparent density inversion at layer ”3” in Figure 67 is not 
significant. 

The net movements of neutrally buoyant particles in the streamers 
and in the invervening layers were observed for several runs. The typical 
net motion (away from the slope in the streamers and toward the slope in 
the intervening layers) is indicated diagrammatically in Figure 68. Aver- 
ages of the net velocity measurements in the streamers and layers are 
listed in the table in Figure 68. Each individual determination of net mo- 
tion that went into the computation of this average involved tracking particles 
for a minimum of 10 wave periods. The differences among the various num- 
bers are not significant; however, a magnitude of about 0. 1 mm/sec in each 
zone is reliable. 



Case 


Average net 
velocities in mm/sec 


Streamer 


Intervening 

Layer 


Critical 
(y = c) 


0. 14 


0.21 


Super 

Critical 

y > c 


0. 17 


0. 13 




Fig. 68. Average net velocities in streamers and intervening layers. 
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(2) Boundary- layer description. Owing to the complicated near- 
bottom motion caused by waves of frequencies at or near critical, it was 
difficult to choose a level that accurately described the physical extent of 
the boundary layer. The height above the bottom of the inflection point in 
the streak curves was selected as an approximate measure of the boundary- 
layer thickness. Admittedly, the choice is somewhat arbitrary; however, 
the change in curvature in the streaks was generally well defined at this 
point during at least the total time of sampling. In any event, this estimate 
of 6 m should be regarded as only approximate. In addition to the values of 
6 , other estimates of the boundary -layer thickness are presented in 

Figure 69 for comparison. * 

The values of fig anc * <5g represent two estimates of the boundary- 
layer thickness that result from scaling considerations. Wunsch (1969) 
found that for co=*co there is a balance in the governing boundary-layer 

C -1/3 

equation over a distance proportional to R -/ . He also showed that the 

equation is nonseparable in the spatial coordinates. This type of balance 
-1/3 

(R ) is represented by 6 0 in column 8. The value of <5« represents a 

* - 1/2 6 

boundary -layer thickness of Stokes type (<5 “R ) and is presented for 
comparison in column 9. Firm conclusions based on quantitative com- 
parisons between the scaling estimates and 6 m are not justified here since 
<5 m is only approximate for each slope station. 

Similar results for a representative run of the supercritical cases 
are shown in Figure 70. For these measurements, the boundary-layer 
motion at each station was very regular, and the dye-streak profiles gen- 
erally resembled those described for the subcritical case. There was no 
evidence of vortex activity with the associated generation of thin streamers 



The estimates of <5j and differ approximately by a factor of 2; this is 
obvious from their definitions. 
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1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


Station 

No. 


X 

c 

(cm) 


h 

(cm) 


r l 

V10- 3 


Rw 


r b 


*i 

(cm) 


A 2 

(cm) 


*3 

(cm) 


5 m 

(cm) 


a 

m 

(cm) 


€ 


1 


15.6 


9.0 


7.86 


68.98 


21. 14 


0. 201 


0.45 


0. 10 


0.362 


1. 18 


0. 12 


2 


24.2 


14.0 


19.03 


92.91 


26. 51 


0. 201 


0. 52 


0. 10 


0.391 


1.37 


0.09 


3 


43.3 


25.0 


60.63 


102. 64 


29.08 


0.201 


0.64 


0. 10 


0.408 


1.44 


0.05 


4 


53.7 


31.0 


93. 2 


25.66 


14.65 


0. 201 


0.68 


0. 10 


0.411 


0. 72 


0.04 


5 


110.0 


33. 5 


109. 2 


7. 15 


5.27 


0.201 


0.70 


0. 10 


0.283 


0.38 


0.02 



*L 

*W 

R 

B 



Nh 2 /V 



o 2 / 

a c o v 

m ' 



(2„/u) 1/2 

hR ,- 1 / 3 



= local wedge Reynolds number 
= local wave Reynolds number 
= local boundary layer Reynolds number 



= (v/N) 



1/2 



a 



FIG. 69. Estimates of boundary layer thickness 
and Reynolds numbers, y = c. 



142 
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4 


5 
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10 


11 


12 


Station 

No. 


X 

c 

cm 


h 

cm 


r l 

V10- 3 


«w 


r b 


6 i 

cm 


6 2 

cm 


* 

cm 


6 

m 

cm 


a 

m 

cm 


€ 


1 


12. 1 


7.0 


4.75 


15. 1 


11.4 


0.22 


0.42 


0. 10 


0.46 


0.61 


0.06 


2 


22.5 


13.0 


16.39 


61.3 


22.4 


0.22 


0.51 


0. 10 


0.45 


1. 23 


0.07 


3 


38, 1 


22.0 


46.95 


30.7 


17.3 


0.22 


0.61 


0. 10 


0.49 


0.87 


0.03 


4 


48.5 


28.0 


76.05 


31.4 


11.8 


0. 22 


0.66 


0. 10 


0.33 


0.88 


0.02 


5 


100.0 


33. 5 


108.86 


1.3 


1. 5 


0. 22 


0.70 


0. 10 


0.21 


0. 18 


0.004 




c 



Nh 2 /y = local wedge Reynolds number 

a 2 oo/y - local wave Reynolds number 
m 

a W 6 /y = local boundary layer Reynolds number 
m m 

(2y/a)) 1/2 

h R l -V3 

h R l _1 ^ 2 = (v/N) 1 / 2 

a k sin a 
m m 



FIG. 70. Estimates of boundary layer thickness and 
Reynolds numbers, y > c. 
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or breaking along the entire slope. Unlike the subcritical ease, with slopes 
of 30 degrees, the displacements of dye streaks and neutrally buoyant par- 
ticles at station 4 did not show any irregularities or anomalous behavior. 
The measurements of the boudary-layer thicknesses at each station are 
considered to be as reliable as those given for the subcritical case 
(Figure 53). 

(3) Velocity and Shear Stress. One-second velocity averages were 
again computed from the streak displacements; for the critical case these 
values are shown only for one vertical level: that described above for the 
measurements of 6 . This level was chosen since it normally was the 

vertical position closest to the slope for which the streak curves remained 
coherent over the entire wave cycle during sampling. The dye streaks 
usually became quite distorted above this level (probably due to the motion 
in the wisps and streamers), and generally the streaks were very faint 
adjacent to the slope (probably due to mixing). These velocity measure- 
ments for the various stations on the slope are presented in Figure 71; 
the experimental conditions at the time these values were determined are 
listed in the upper right-hand section of Figure 69. 

The values of velocity in Figure 71 are plotted along the vertical 
axis in cm/sec, and the successive seconds of time during the wave cycle 
are shown along the horizontal. The phase has been adjusted artificially 
so at time t = 1 the velocity is nearly zero. The average distance above 
the bottom for which these measurements were taken is indicated for each 
station. This distance is approximately equal to S m (Figure 69). A sum- 
mary of the maximum and RMS values of velocity is given in Figure 72 for 
each station (including the velocity variations on the flat bottom not shown 
in Figure 71). 
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The last column of Figure 72 gives rough estimates of the maximum 
bed shear stress, 



I r I max 
1 o 1 






1 ul max 

6 



where ~jl = 0. 01 gm/cm-sec, that were computed at the level of the inflec- 
tion point (if = <5 m ) described above. This kind of estimate assumes a 
laminar velocity distribution that remains essentially linear with distance 
above the bottom over the wave cycle. This assumption is not accurate 
for two reasons: (1) the actual profiles of velocity are probably not linear, 
and (2) the observed dissipation of the dye streaks below 5 m indicates pos- 
sible turbulence in this region. The estimates of It I are intended 
only as approximate lower bounds on the local maximum bed shear stress; 
the actual values are probably larger due to a nonlinear velocity distribu- 
tion with superimposed turbulent fluctuations near the bed. 

Fig. 72. Summary of velocity and bed shear stress 

values for critical case (y^ c) measurements. 



Station 
(see Fig. 11) 


1 u 1 max 
cm/sec 


U RMS 

cm/sec 


l 7- o Imax 
dynes/cm^ 
(xl0+2) 


cm 

<5m 


1 


0.65 


0.38 


1.80 


0.36 


2 


0.68 


0.44 


1.74 


0.39 


3 


0.72 


0.45 


1.76 


0.41 


4 


0.35 


0.25 


0.85 


0.41 


5 


0. 18 


0. 13 


0.64 


0.28 



The velocity measurements show that the magnitude of the maximum 
oscillatory motion at the selected level of measurement is nearly uniform at 
the middle and upper slope stations. The value of | u | max on the lower slope 
(station 4) is smaller than that farther upslope (stations 2, 3, and 4) by about 
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a factor of 2. The increase in maximum velocity between the flat bottom 
and station 1 is approximately a factor of 4. This is considerably lower 
than that found between similar locations for the subcritical-case measure- 
ments shown in Figure 59. The measured bed shear stress (as defined 
above) at station 1 is greater by only a factor of 3 than that for the flat 
bottom. A striking feature in the results is the smooth sinusoidal appear- 
ance of the velocity- time curves in Figure 71. The apparent lack of irreg- 
ularities in these curves suggests that in spite of the vortex-streamer 
production, the basic harmonic motion of the input waves is still quite 
recognizable. 

The Lagrangian velocity measurements for the supercritical case 

are represented by the results in Figures 73a and 73b. These velocity- 

depth profiles are analogous to those shown in Figures 54 through 58 for 

y < c. Only the results at two stations (1 and 4) are shown here to illustrate 

the regularity of the motion. The absence of the vortices and thin streamers 

for this frequency (u> < u> ) was also accompanied by a return to an apparent 

c 

periodic motion of the dye streaks and particles near the bottom. No wave 
breaking was observed near the corner. The maximum and RMS velocity 
values that were taken from selected levels in the velocity-depth curves at 
each station are listed in Figure 74. The maximum value of bed shear 
stress is also shown in the last column of this figure. 

The velocity curves at all stations are quite symmetrical and do not 
indicate any measureable net motion; the measured net displacement of 
neutrally buoyant particles over several wave periods corroborates this 
finding. It is also evident that 1 u I at the selected levels increases 
only slightly from station 4 to station 2, and then decreases at station 1, 
which is farthest upslope. This behavior contrasts sharply with the large 
increase in | u l max station 1 for the high-frequency waves (y <c, case A). 
The estimated bed shear stress at the slope locations increases by a factor 
of 4 or 5 from that over the flat bottom. In this case, <5 in Equation (4-7a) 
was chosen as the smallest jt for which the dye streak data were considered 
to be reliable (see first value of rf for each station in Figure 74). 
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FIG 73b. Velocity - depth profiles at station 
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FIG. 74. Velocity and bed shear stress at selected levels 



|t 0 I max 



Station 

No. 


V 

cm 


| u | 

1 'max 
cm/sec 


1 u Irms 

cm/sec 


(Xio 2 ) 

dynes/ c m 


1 


0. 11 


0. 18 


0.08 


1.64 




0.26 


0.26 


0. 14 






0.38 


0.29 


0. 17 






0.45 


0.28 


0. 17 




2 


0. 16 


0.32 


0.22 


2.00 




0.30 


0.44 


0.30 






0.45 


0.48 


0.33 






0.66 


0.50 


0.33 




3 


0. 13 


0.20 


0. 13 


1. 54 




0.24 


0.30 


0.20 






0.38 


0.35 


0.24 






0.53 


0.36 


0.24 




4 


0.13 


0.25 


0.17 


1.92 




0.24 


0.35 


0.23 






0.37 


0.38 


0.25 






0.51 


0.35 


0.23 




5 


0.12 


0.05 


0.03 


0.42 




0.21 


0.07 


0.04 






0.32 


0. 11 


0.06 






0.46 


0. 11 


0.07 
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DISCUSSION 

The experimental results can be summarized briefly. The measure- 
ments of the wave motion in the subcritical range y <c indicate that the 
linear theoretical solutions (Wunsch, 1969) accurately represent the data 
for the higher frequency waves in the interior and in the boundary layer 
over the slope for e less than approximately 0.5. The data for the lower 
frequency waves in the subcritical range do not agree with these solutions 
as well as the data for the higher frequency waves. For both of these wave 
types, no significant reflections were measured; the wave energy dissipates 
primarily in a narrow upslope zone. The higher frequency waves break 
turbulently in this zone, where mixing generates laminae of various thick- 
nesses that protrude back into the interior. By contrast, the lower fre- 
quency waves form surges that dissipate in runup near the corner. Large 
intensification in the velocities parallel to the slope were measured for all 
waves in this range. 

The critical case, y ^ c, is characterized by less violent breaking 
near the corner. Regularly spaced line vortices form along the slope dur- 
ing the upslope motion of the water particles. The mixing in these vortices 
generates thin streamers of fluid that penetrate horizontally back into the 
interior over most of the depth. Intensification of the near-bottom oscillatory 
velocities was also evident; however, the flow pattern in the vortices com- 
plicated the kinematics near the bottom. 

In the supercritical case, y > c, no breaking was seen near the 
corner. The vortices and streamers persisted until the wave frequency 
reached low values, below which the vortex production subsided. The 
motion in the boundary layer was very regular in the absence of the insta- 
bilities. In the fluid interior the measurements show that the wave motion 
intensifies near the critical characteristics, whose position is given by 
z = -cx. 
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5. SEDIMENT MOVEMENT BY INTERNAL WAVES 
GENERAL DISCUSSION 

This chapter attempts to assess the role of shoaling internal gravity 
waves in the initiation of bottom-sediment movement in the ocean, in the 
light of the theoretical and experimental results discussed in the previous 
chapters. Although the concept of internal waves as an effective geological 
agent in the ocean is not new, very few published studies have dealt speci- 
fically with the erosion, transportation, or deposition of sediments on sub- 
marine slopes. Chapter 1 has provided some historical background on this 
subject. Several explanations for the apparent lack of attention are possible: 
(1) our knowledge of the mechanics of sediment movement beneath the more 
familiar oceanic surface waves is still incomplete; (2) the wide spatial 
separation of available data on bottom sediment types on the continental 
margins hampers understanding of local sedimentary processes; (3) reli- 
able measurements of oceanic internal waves, particularly those taken in 
conjunction with a sediment-sampling program, are few; (4) until very 

recently, theoretical and experimental studies of internal waves passing 

* 

over a bottom with variable depth were virtually nonexistent. In spite of 
these problems, there have been inferences to the potential importance of 
this process (Lafond, 1962; Hulsemann, 1968). There exists a body of 
internal-wave measurements (see Chapter 1) which indicate that internal 
waves propagate shoreward normal to the bottom contours. (Wunsch, 1969, 
has also shown theoretically that these waves can refract, in the sense that 
their crests tend to align parallel to the bottom contours.) These measure- 
ments, together with the results of the previous chapters, promote the fol- 
lowing discussion about the interaction of these waves with bottom sediment. 



% 

Defant (1961) gives two illustrations of the experimental work of Zeilon 
(done in 1934) in which interfacial waves incident onto a model shelf-slope 
configuration were studied. 
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The first section in the following discussion presents various 
examples of internal -wave measurements in the ocean and discusses quali- 
tatively some features of sediment distribution on the continental margins, 
particularly off the northeastern United States. The purpose of this section 
is to provide some background material on ocean conditions relevant to the 
later discussion and to demonstrate the wide range of values of the various 
hydrodynamic and sediment properties. The second section develops a 
criterion that describes the initiation of bottom -sediment motion. This 
criterion is basically an extension of the analysis developed by Eagleson, 
et ad (1958) for incipient motion of discrete bottom-sediment particles in- 
duced by surface waves. The analysis for shoaling internal waves over a 
constant bottom slope is then specialized for the case of small slopes 
(y <<c), which is shown to be particularly appropriate for high-frequency 
internal waves in the ocean. Another criterion (the Shields curve, with a 
later extension by Vanoni, 1964), which represents the more traditional 
approach to the problem of incipient motion of bottom sediment and is based 
primarily on dimensional analysis and empirical results, is also discussed. 
This latter criterion is shown to be approximately equivalent to the first. 

The analysis is then applied to specific oceanic measurements, and 
the results are discussed in a final section, in light of the various limita- 
tions. Qualitative remarks concerning the effect on sediment movement by 
breaking internal waves and critical slope conditions (y =* c) are also given 
in the last section. Throughout the subsequent discussion, consideration of 
the interaction between internal waves and sediment particles is only an 
initial attempt to evaluate the possibilities of incipient motion induced by 
this interaction; the actual importance, or even existence, of this process 
under real conditions must await future field studies. 

Oceanic Conditions 

Chapter 1 provides several references to actual measurements of 
internal wave motion in the ocean, including studies of the directional 
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properties of these waves (for example, Summers and Emery, 1963). In 

most cases the measurements were made either by recording fluctuations 

of temperature or horizontal velocity at fixed positions from a moored buoy 

or anchored ship or by sampling the vertical distribution of temperature 

with devices like the bathythermograph. Since the frequency content of the 

internal wave field is not known a priori, a common means of distinguishing 

the contributions to the measured changes by the various frequency compo- 

* 

nents is the energy spectrum. Figure 75 illustrates kinetic energy spectra 
at several levels computed from horizontal velocity data taken at site "D’ 1 
(location shown in Figure 76) and described by Fofonoff (1969). The large 
values of spectral energy density at the inertial and tidal frequencies are 
typical of this location. Frequencies and amplitudes of vertical tempera- 
ture fluctuations have been measured by others (Chapter 1). Figure 77 
demonstrates the wide range of amplitudes and frequencies represented in 
previous measurements of internal waves in the ocean. The measurements 
by Gaul (1961) and Ufford (1947) taken over the outer continental shelf and 
continental slope suggest that wave amplitudes of 1 meter to 5 meters are 
not unreasonable for the higher frequency internal waves. 



Time series of measurements of vertical velocity or temperature fluctua- 
tions can be used to compute frequency spectra of potential energy. Mea- 
surements of horizontal velocities like those taken at side "D" (location 
given in the following text) are the basis for kinetic energy spectra. 

Fofonoff (1969) discussed the spectral ratio of potential energy to kinetic 
energy in terms of a linear internal wave model and found that this ratio is 
approximately unity for low frequencies (near inertial frequency), but ap- 
proaches infinity for co -*■ N. He noted that this latter behavior is probably 
caused by the exclusion of the advective terms in the linear model. Voor- 
his (1968) measured vertical internal motions near site "D" with neutrally 
buoyant floats. After comparing the frequency spectra of potential energy 
and vertical kinetic energy (obtained from the measurements) with the spec- 
tra of horizontal kinetic energy computed from the site "D" records, 

Voorhis found that there is an approximate equipartition of potential and 
kinetic energies for frequencies between inertial and Brunt-Vaisala. He 
interpreted this equipartitioning to indicate the predominance of internal 
wave motion in the measurements. 



KINETIC ENERGY cm /sec / cph 
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Fig. 75. Frequency spectra of horizontal kinetic energy 
density at site "D". 
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Bathymetry off southern New England. 
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FIGURE 77 

INTERNAL WAVE MEASUREMENTS ON CONTINENTAL MARGINS 



Source 


Water Depth 
(meters) 


Wave Height 
(meters) 


Wave Period 


Celerity* 

(knots) 


1 . 


Lafond 

(1963) 


20 


7 (max.) 


7 min. to 
12 hr. 


0. 11-0.6 


2. 


Boston 

(1964) 


20, 33 


8 


12 hr. 




3. 


Gaul 

(1961) 


62 


1 - 3 


10 min. 


0.5 - 1.5 


4. 


Ufford 

(1947) 


40, 100 


1 - 30 


9 min. - 
2 hr. 


0.08-0.68 


5. 


Emery 

(1956) 


1350 


130 - 200 


12 hr. 


(standing) 


6. 


Lee 

(1961) 


20 


3 - 7 


5-15 min. 


— 


7. 


Summers and 

Emery 

(1963) 


10 - 1500 


30 


12 hr. 


1 (shallow) 
7 (deep) 



*In all cases, the dominant direction of propagation was shoreward. 
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In order to describe internal wave motions in the ocean, an estimate 
of the vertical distribution of Brunt-Vaisala frequency is needed. The aver- 
age vertical distribution of Brunt-Vaisala frequency for locations seaward 
of the continental slope (near site ”D", for example) is typically character- 
ized by relatively large values of N in the seasonal thermocline, if present, 
and in the deeper main thermocline. Figure 78 shows measured values of 
Brunt-Vaisala frequency N versus depth for stations near site ”D”; the hydro- 
graphic data were provided by Volkmann (personal communication). An 
idealized distribution of Brunt-Vaisala frequency with depth is also shown in 
this figure. The idealized case consists of three distinct vertical layers; in 
each layer the density gradient is assumed to be constant (i. e., N ^ constant). 
This hypothetical vertical distribution of N is later used with the idealized 
topography also shown in Figure 78 as a basis for various models of sedi- 
ment movement induced by internal waves. 

The generalized topographic cross-section in Figure 78 is not con- 
strued to be representative of any particular area, but is taken as an 
idealized case in which the slopes are linear and the boundaries between 
the three provinces (shelf, slope and rise) are well defined. Each idealized 
spatial province in this figure is similar to the experimental conditions 
described in previous chapters in the sense that N and y are constant in 
each. 

The actual bottom topography of the various provinces comprising 
the continental margins of the oceans not only varies markedly with geo- 
graphic position but also can be quite complicated locally. Featureless 
and even bottom slopes are probably rare. In his summary of the continental 
shelves and slopes of the world, Shepard (1963) remarks that it is exceptional 
to find a shelf that has a constant slope out to the shelf edge; constancy of 
slope usually implies widely spaced sample points. The values for bottom 
slopes shown in Figure 78 are recommended averages taken from Shepard 
(1963) and Heezen, et al (1959). More detailed descriptions of the physiog- 
raphy, structure, and sedimentology of the continental margin off the east 



shelf 
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FIG. 78. Idealized cross section of continental nargin. Actual and idealized distribution 
of Brunt-Vaisala frequency are also shown. 
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coast of the United States is presented by Uchupi (1963, 1968), Emery (1966) 
Pratt (1968), Schlee (unpublished manuscript), and Emery, etal(1969). 

At this point it is instructive to compare average values of the 
oceanic slopes y with hypothetical estimates of the slopes of the internal 
wave characteristics c in order to determine which case (y <c, yac, or 
y >c) in Figure 1, if any, is predominant. Figure 79 summarizes the 
relative values of c and y for typical values of Brunt-Vaisala frequency N, 
wave frequency co, and a constant value of Coriolis frequency (f = 0.042 cph). 
The results show that the continental shelves are approximately critical 
(i. e., y a c) for internal wave frequencies near inertial for each of the values 
of N shown; for frequencies greater than the frequency of the semi-diurnal 
tide o> T , the shelf is subcritical. 

The assumption that the average gradient of the continental shelf is 
small compared to the slope of characteristics c for normally incident 
internal waves (i.e., y <<c) appears to be reasonable for high-frequency 
internal waves. For example, if w>0. 1 cph and N > 1. 0 cph, or if o>>0. 5 
cph and N = 6. 0 cph, then y <<c is approximately valid according to 
Figure 79. The average gradient of the continental slope might be critical 
for a relatively wide range of internal wave frequencies depending on the 
local values of N and f. For example, in Figure 79 for N in the range 
considered (0. 5 cph <N <6.0 cph), the range of critical frequencies is 
approximately 0. 045 cph < co <0. 03 cph; this range includes the semi- 
diurnal tide. The average gradient of the continental rise is approximately 
critical for frequencies between inertial and tidal (semi-diurnal) for 
0. 5 cph < N < 3. 0 cph. 

It is tempting to speculate about the consequences of y > c for 
oceanic conditions; this is done later and only touched upon here. The ex- 
perimental results show that for internal wave frequencies equal to or 
slightly less than critical, vortex instabilities form along the slope. If 
this occurs along submarine slopes, the mixing in these vortices, with 
consequent formation of thin streamers that might advect fluid horizontally 
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FIG. 79. Characteristic slope c (ordinate) computed for various 
frequencies t-*j- and selected Brunt Vaisala frequencies /V 
Average oceanic slopes Y are also shown. 
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back into the interior, might provide a means of stirring the fluid near the 
bottom and moving the very fine suspended matter seaward. This agitation 
might inhibit or prevent deposition of fine suspended matter, and thus main- 
tain a relatively turbid and well-mixed fluid layer near the bottom. How- 
ever, it is dangerous to extrapolate the laboratory observations to unobserved 
natural phenomena. 

The intention of the above discussion of average gradients of the 
oceanic slopes was to present the gross features of these slopes in relation 
to an idealized hydrodynamic structure (Figure 79); local bottom features 
such as animal markings, sediment ripples, depressions, and rock outcrops 
were intentionally left out of consideration. It is realized that these additional 
features might locally affect the internal wave motion by causing complicated 
reflection patterns of the internal wave characteristics (Longuet-Higgins, 

1969) and by influencing the flow in the boundary layer. 

The actual scales of bottom roughness on oceanic slopes are quite 
variable. Spatial variations of these scales have been illustrated by several 
excellent photographic and echo-sounding studies of submarine topography. 

A series of camera stations spaced two miles apart on a line extending south 
of Martha's Vineyard and lying approximately normal to the bottom contours 
was described by Northrup (1951). Sediment samples were taken at most of 
the camera stations. Northrup noted that there was good evidence for the 
action of bottom currents on the outer shelf. Others (for example, Uchupi, 
1963) have described topographic features such as sand waves with heights 
of 5 to 25 meters on the continental shelf. Owen and Emery (1967) cited 
photographic evidence for intensive erosion of the continental slope south of 
of Martha's Vineyard by strong bottom currents. They found several un- 
usual features in their photographs including small masses of sediment that 
protruded a centimeter or so above the general level of the bottom; in this 
instance they described the overall appearance as "sand-blasted". Bowin, 
et al (1967) made an extensive photographic and sediment survey of the 
western and southern slopes of Plantagenet Bank (near Bermuda) in water 



162 



depths of 60 meters to 2000 meters. The height scales of the various sand 
ripples that are clearly shown in some of their stereo -photographs appear 
to be on the order of 10 cm. (The ripples were observed at various depths, 
but especially between 950 and 1200 meters.) 

The effects these local variations in bottom roughness actually have 
on the internal wave motion have not been determined. In fact, what role 
internal waves might have in forming or maintaining features such as ripples 
or sand waves is uncertain, and must await field studies of this problem. 

The analysis and discussion in the subsequent sections deals only with the 
simpler problem of a linearly sloping bottom, mainly for reasons of analy- 
tical tractability. It is probably true that the wide variations in the hydro- 
dynamic and topographic features of the ocean include this case; its consider- 
ation might provide some insight into more complex situations. 

A general discussion of the distribution of sediments on the various 
continental shelves and slopes is given by Shepard (1963). For the purposes 
of this discussion, only a brief summary of the sediment types in the area 
shown in Figure 76 will be given. Uchupi (1963) and Schlee (unpublished 
manuscript) presented detailed descriptions of the bottom sediments in 
this region. In his report Schlee noted that although the sediments on the 
shelf are dominantly quartzose sands that are moderately sorted, many of 
the sediment samples from the continental shelf (about 20 percent) showed 
sizable silt-clay fractions. Schlee also noted that there is some evidence 
for a seaward reversal in the median sand sizes (i. e., fine to coarse in a 
seaward direction) at limited sections of the shelf break southeast of Cape 
Cod. This seaward reversal has been found by others. 

Several studies have indicated that the sediments on the outer shelf 
shown in Figure 76 are predominantly Tertiary, and that modern sediments 
are either lacking or the result of reworking of the older material (Donahue, 
et al 1966). The dynamics of the by-passing of the shelf by much of the silt 
and clay of modern sediments, particularly on that part of the shelf below 
surface wave base, are not well known. Emery (1966) states that waves 
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and currents on the continental shelf produce enough motion to cause much 
of the silt and clay of modern sediments to travel across the outer shelf 
and come to rest in deeper waters. It is suggested later in this chapter 
that near-bottom velocities induced by internal waves might be partly re- 
sponsible for the sediment transport in this region. 

Pratt (1968) described the continental slope off the eastern coast of the 

United States as a complex feature whose surface is generally more irregular on 
the upper parts . In the area under consideration the median diameter of the sur - 
face sediments is in the silt range; quartzose sands are locally predominant on the 
upper slope (Schlee, unpublished manuscript). The continental rise has been de- 
scribed as a depositional apron sloping seaward from depths of about 2000 m to 
about 5000m at a gradient of between 1:100 and 1:700 (Heezen, et al 1959). The 
surface is covered with fine well-sorted sands and coarse silts interbeddedwith 
clays. Emery (1966) suggests that these sediments were deposited by turbidity 
currents; however, Heezen, etal (1966) maintain that ocean bottom currents 
flowing parallel to the contours are responsible for most, if not all, continental 
rise sedimentation. 

Actual sediment data that are later compared with the analytical 
results were obtained from the work of Hollister (personal communication), 
Stetson (1939), Pratt (1968), and Emery and Ross (1968). The geographic 
locations of the sampling sites are shown in Figure 76. These results are 
particularly useful, since the sample locations are on lines extending off- 
shore near site "D". The hydrodynamic data at and near site "D" provide 
convenient inputs to the models considered in a later section in which 
incipient-motion criteria are discussed for the various sediment data. 

Before proceeding to the analysis of sediment movement by internal 
waves, some consideration of the boundary layer along oceanic slopes would 
be worthwhile. Studies of sediment movement in rivers and beneath surface 
waves have indicated the importance of the flow characteristics (laminar or 
turbulent) in this layer in determining the motion of bottom sediment. Linear 
solutions due to Wunsch (1969) for the boundary layer beneath internal waves 
on a sloping bottom were discussed in Chapter 4. Wunsch (1970) also 
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considered theoretically oceanic boundary mixing for quasi-steady conditions 
in a stratified fluid adjacent to a sloping boundary, and showed that a mean 
vertical velocity is induced at the sloping boundary to satisfy the no-flux 
condition. He found that the effect is confined to a boundary layer of thick- 
ness 

6 = (vk) 1/4 /N 1/2 

where 

v is the eddy momentum coefficient 

k is the eddy heat coefficient 

N is the Brunt-Vaisala frequency, as before. 

He also pointed out that the various estimates of the mixing coefficients 

2 

v , k in the ocean are unreliable, and showed that for v , k ~ 1 cm /sec, 

-3 -1 

and N = 2 X 10 sec , then 6 ~ 20 cm. 

Attempts to define v for various oceanic processes are abundant 

(Defant, 1961); however, near a sediment -strewn bottom that is disturbed 

by periodic and quasi-steady currents, such a definition is elusive. Based 

on measurements of suspended matter near the bottom of the continental 

slope east of Chesapeake Bay (Ewing and Thorndike, 1965), Ichiye (1966) 

computed an eddy diffusivity from a linearized version of the Fickian dif- 

2 

fusion equation. He arrived at a figure of 0. 12 cm /sec for the eddy 
viscosity. Others have shown photographically that the bottom along the 
submarine slopes can be quite turbid; high concentrations of suspended mat- 
ter in this near-bottom zone might increase both the viscosity and the den- 
sity of the fluid. By contrast, many photographs show remarkable clarity 
and suggest low concentrations of suspended material near the bottom. It 
therefore appears that the local conditions along oceanic boundaries are 
sufficiently variable to preclude a single, all-inclusive value of v • The 
above discussion suggests that many processes are operative along the 
oceanic slopes (currents, animal stirrings, etc.); internal waves are only 
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one possible kind. The actual viscosity near the bottom must be the result 
of the composite effects of the processes. 

INCIPIENT MOTION CRITERIA 

The various approaches to the formulation of analytical criteria 
that adequately describe the initiation of motion of bottom sediment due to 
water waves are limited owing to the complexity of the process. Several 
authors (Inman, 1963; Raudkivi, 1967; Abou-Seida, 1965) have summarized 
these approaches and discussed some of the limitations of each method. 

Prior to the development of an analytical model, it might be helpful 

to present a dimensional analysis of the problem of initiation of sediment 

motion induced by internal waves, in order to show the large number of 

independent variables and dimensionless groupings of these variables that 

are possible. The specific problem is to relate functionally the incipient 
* 

diameter of the movable sediment particles to the pertinent hydrodynamic, 
geometric, and sediment variables. It is assumed that internal gravity 
waves of a given mode number in a continuously stratified fluid propagate 
over a bottom slope that has a sedimentary bed. The density stratification 
is approximately linear and stable, and the bottom slope is constant. The 
bed is assumed to consist of noncohesive particles with mean diameter D^; 
the roughness is arbitrarily described by a length scale < . Figure 80 
shows the idealized geometry and particle relationships; this figure is 
similar to that given by Eagleson and Dean (1959). 



^Incipient diameter is the largest size within a local size-frequency distri- 
bution of particles that can be placed into motion by the action of the hydro- 
dynamic and gravitational forces. The condition of incipient motion induced 
by surface waves on beaches has been described as follows. At some point 
on the offshore slope, the instantaneous hydrodynamic forces may become 
sufficiently large to cause static instability of bed-sediment particles of a 
given size. The diameter of this size is referred to as the incipient diam- 
eter (Johnson and Eagleson, 1966). 
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symbol 



Definition sketch of sediment particles; forces and 
geometry. 



name 



definition 



point of 
application 



VM 



submerged 

weight 

fluid 

resistance 



pressure 

force 



virtual mass 
force 



ttD /6 p f g (s g -s f ) 

gravity - buoyancy 

dynamic pressure and 
viscous shear acting 
over surface of particle 

7 rD ^/6 p- a 

f — o 

instantaneous pressure 
gradient force present in 
fluid whether or not 
particle is present. Equal 
to the inertial force of 
the displaced fluid. 

C M *° 3 / 6 Pf -f 

force required to accel- 
erate fluid due to 
presence of particle 



center of 
mass 



top 



center of 
mass 



top 
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Eleven independent variables are selected a priori as being important 
in determining the incipient diameter D. for this problem. The functional re- 
lationship can then be expressed as 



(The symbols are defined in the front section of this paper.) By dimensional 
analysis the functional relationship can be expressed equivalently in terms of 
eight independent dimensionless groups; 



where s = s /s.. 

s' f 

The first four expressions on the right in Equation (5-2) relate to the 



local wave field ; the remaining groups relate either to the wave-sediment 
interaction or to the sediment properties alone. Dimensionless groups (5) 
through (8) resemble those developed by Carstens, et al (1969) in his analy- 
sis of incipient motion due to surface waves. The large number of dimension- 
less groups indicates the complexity of the process and the potential difficulties 
in developing accurate dynamic models of this process. Several other vari- 
ables have purposefully been omitted; these include measures of sorting and 
packing. These factors would increase the complexity of the process. 

If the problem is simplified to consideration of the balance of forces 
acting on a single, protruding bed particle such as shown in Figure 80, an 



a/h is proportional to the local Stokes parameter e, since k ^ cn 7 r/k anc * 
e = ak (Chapter 4). 



D i - f(a, u), N, n, p f , y, h, y, D m , p g , k) 



(5-1) 



D./D 
r m 



f/a/h, co/N, n, y/c, acoD^/y , a co/[ (s - 1) gDj 1//2 

V) (2) (3) (4) (5™ (6) 



(5-2) 
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expression can be derived for the incipient diameter in terms of the local, 
wave-induced velocity field in the viscous boundary layer, certain resistance 
coefficients, bottom slope, and sediment properties. Furthermore it is 
postulated that incipient conditions hold when the sum of all hydrodynamic 
and gravitational moments acting on the particle at any instant is equal to 
zero. Such a condition can be represented (similar to Eagleson and Dean, 
1959) by 

V M = 0 = (F f sin 0 sin <f> ± F p cos 0(1 + cos 4) 

± Fy M (1 + cos <f>) ± Fp cos (f) - F g sin (<f> +a)) ^ 

Brief definitions of each force and its assumed point of application 
for this model are given in Figure 80. The points of application follow from 
the discussion in Ippen and Eagleson (1955) and Eagleson and Dean 0959). 
Ippen and Eagleson (1955) give a more detailed discussion of the various 
forces and their representation. The fluid resistance can be rewritten in 
terms of components normal to and parallel to the bed; these components 
are commonly referred to as lift and drag, respectively (Schlichting, 1955): 

F l = F f sin 0 = C L p f ttD 2 /8 u f 2 

F d = F f cos 0 = C D p f 7 T D 2 /8 u f 2 

To obtain an expression for the incipient diameter from Equation 
(5-3), the following assumptions are made. 

(1) The magnitude of the acceleration is small compared to the 
magnitude of the effective velocity (|aj, | << | u f | ) . If the advective 
accelerations are neglected (this is a linearized model), then the instanta- 
neous accelerations are approximately equal to co times the instantaneous 

velocity. Their neglect is reasonable for internal waves in the ocean (say 

m-3 -U 
co ^ 10 sec ). 
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(2) The lift coefficient is approximately equal to the drag coefficient. 
The actual relationship used in this model (C T » 0. 85 C^) was obtained from 
Chepil (1958). It is recognized that this relationship was found for particles 
in a turbulent boundary layer; its applicability to laminar conditions has not 
been verified. 



Applying assumptions (1) and (2) above and the lift and drag expres- 
sions for the fluid resistance in Equation (5-3) and solving for D, the incipi- 
ent diameter is expressible as 





g 




(5-4) 



where 



± (1 + cos <ft) + 0. 85 sin 
( s s/ s f - 1) sin Wra) 



The sign convention necessitates the use of absolute value for r; if u f is 
directed upslope the choice of the bottom signs would give a negative value 
of r. The value of represents the largest diameter that can be moved 
for particular sediment and internal-wave conditions. 

The quantity r is dependent on<£, s /s f , and a. Since emphasis is 
placed here on noncohesive particles in sea water, the choice of s g = 2. 65 
(quartz) and s^ = 1. 027 appears reasonable. The density variation in the 
water column over the continental margins off the northeastern United 
States is < 1 percent. The angle of repose <£ for both well-rounded and 
very angular grains of sizes >0.3 mm was shown graphically by Albertson, 
et al (1963). There is strong evidence from the nature of the curves for 
rounded and angular grains of small diameters (<0.3 mm) that <£ approaches 
30 degrees (± 3 degrees). Eagleson, et al (1958) showed that for the idealized 
situation similar to Figure 80 this angle can be derived geometrically: 

0.866 

[(D/k) 2 + DA - 1/3] 1/2 



tan <£ = 



(5-5) 



170 



If the diameters of the bed grains are approximately equal to the equivalent 

diameters of the roughness (D/k * 1), tan <f> = 0. 53 or $ =* 28 degrees. If 

these values of s , s-, and <j> are used and a is assumed to be small, such 
s f 

that a << <f>, then | r ] =* 2. 01 for upslope motion and | r | =* 3. 08 for down- 
slope motion. 

Consequently, on small slopes (a << </>)> 
u f 2 C D 

D. = 2.31 (downslope incipient motion) (5-6a) 

^ S 

and 

u f 2 C D 

= 1.51 (upslope incipient motion). (5-6b) 



The problem now remains to obtain a value u^ in Equation (5-6a) or 
(5-6b) that is the maximum velocity which effectively acts at the uppermost 
point of the grain surface (rj" = D.). The boundary-layer velocity at this 
point can be expressed as the real part of Equation (4-24) 



u f = ^ [ { F sin6j - G sine^} - e ^i^ 6 

X {F sin(0 1 + D^/s) - G sin (0 2 + D^s)}] 

..(-HjH-r 

\cu - 1ST sin a/ 

where 

a - 2ni A - 

q CnA ’ c - r 

0 1 = qCn(BJ) + w t , 6 ^ = q£n(Pf ) + cut 



(5-7) 



(5-7a) 



For the case of small slopes in the sense y << c, Wunsch (1969) has shown 
that q « Furthermore, for distances well away from the corner 

(cx > > |z j)| the coefficient A can be written in terms of the input wave 
motion 
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a o> 




also k 

o 



cnjr 

h 



o 



(5-7b) 



where a and k are the amplitude and wave number of the internal waves 

oo 

over the flat bottom, respectively. This expression is true for an input 
area whose top and bottom are horizontal and rigid. If the condition that 
the vertical velocity of the wave motion at the top and bottom of the input 
channel is approximately equal to zero, then Equation (5-7b) is approximately 
true for an input channel with nonrigid upper and lower boundaries. It is 
important to note that these approximations of q and A imply a smooth 
transition in flow conditions from the input region to that over the slope 
(i. e., in two dimensions this implies a smooth matching of stream lines at 
the transition between the two regions). 

For conditions of small slope (y << c) and a smooth bottom transi- 
tion, several simplifications are possible in Equation (5-7). For small 
slope, y « a ~ sina and for high enough wave frequencies, co/N>> sin a. 

In addition, 



M =* c; N » c 




so that 




a Q om7r 



sine g 



-D/s 

e sin (e 3 + D./s ) 



where 

e 3 = q£n(c£) + cot 



6 « (-) 
CO 



1/2 



Since k = -r — ^ and h = -y x «-yi for small y , 
oh ' ' 

o 



then 



a o w 

u f “ ^h/h Q ) 



sine 3 



"D /S 

e sin (0 3 + D./6 ) 



(5-8) 



(5-9) 



with n = 1. 
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After a short amount of manipulation it can be shown that the maximum 

value of the terms inside the brackets on the right side of Equation (5-9) 

can be approximated as \J 2 D./s to order (D./§) . For waves of 1-hour 

-3 1 1 -9 2 

period, co 1.7 X 10 rad/sec, and for v =* 10 cm /sec, 6 =* 10 cm. 
Typical sediment diameters on the outer shelf are on the order of 10 ^ cm, 
indicating that the above approximation is reasonable to about one percent. 
Since for incipient motion conditions it is the maximum value of the velocity 
at rj = that is the important quantity, then from Equation (5-9) 



a co 

o 



U flmax “ c(h/h ) yj2 D i^ 6 



(5-10) 



for the following assumptions: 

(1) small slopes such that y << c and a <<<f>; 

(2) D. << 5 . 

It is also worth recalling here that the results of the laboratory experiments 
described in the earlier chapters for the subcritical case (y < c) suggest that, 
at least for the conditions tested, the velocity field in the boundary layer was 
adequately described by the linear theoretical solutions along the slope to the 
point of breaking. For any particular case in which y <<c is not valid, but 
y < c still holds, the complete expression for u^ (Equation (5-7)) must be 
used in Equation (5-4). However, it has been shown (Figure 79) that y << c 
is reasonable for several average oceanic conditions. 

If Equation (5-10) is applied to Equation (5-6a), then 

2 

g 5 (downslope incipient /'K-m 

motion) ' 

This equation shows that in order to determine some estimate of boundary 
layer thickness 5 (or equivalently, viscosity y) must be made. Owing to the 
difficulties in obtaining oceanic measurements of this quantity, this choice is 
not simple. However, a lower bound can be estimated: the value of 6 obtained 



D. = 0.216 

l 



c (h/h o ) 
. a o" 
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from Equation (5-7a) with v = where is the kinematic viscosity 

(i/j^ « 10 2 cm 2 /sec). The variation in for sea water is shown by 

Higgins (1962); typically v . varies from about 0.008 cm 2 /sec at 30°C to 
2 K 

about 0.015 cm /sec at 0°C (salinity about 35% 0 )- Although coefficients 
of viscosity have been suggested by other authors to describe mixing 
processes that have scales larger than molecular, the current state of 
knowledge of the oceanic boundary layer induced by internal waves does not 
justify an a priori assumption of v >y^. The problem of selecting a repre- 
sentative value of v was discussed at the end of the previous section in this 
chapter. The basis for other choices of viscosity here lies in an "intuitive 
feeling" that natural oceanic boundary layers are turbulent. The actual 
viscosity is probably a result of several processes at work, including in- 
ternal waves. It is worth recalling, however, that the experimental runs 
with linear, smooth slopes showed no appreciable net transport of fluid in 
the boundary layer outside of the breaking zone for y <c. In the discussion 
of net Lagrangian motion in Chapter 4, this observed absence of net trans- 
port was interpreted to mean that there was no measurable mixing induced 
by the internal waves along the slope up to positions near the breaking zone. 
The turbulent dissipation of the high-frequency waves for y <c in the break- 
ing zone and the mixing induced by the vortex instability for y> c are two 
possible instances for which a choice of v >i/^ might be made, at least for 
the laboratory conditions. 

The other apparent difficulty in the direct application of Equation 
(5-6a) or Equation (5-6b) to oceanic conditions is the choice of the drag co- 
efficient Cjy Considerable experimental evidence shows that C^ is largely 
dependent on a suitable Reynolds number for solid bodies immersed in a 
moving fluid (Prandtl, 1952; Batchelor, 1967). Usually the relationship be- 
tween Cp and Rp is given for a steady flow past a smooth and regular body 
such as a sphere in an unbounded fluid. In the case of an irregularly shaped 
particle resting on a rough bed in an oscillatory flow, the relationship be- 
tween C^ and Rp is not well known. Eagleson et al (1958) showed 
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measurements from which the drag coefficients for spheres rolling down a 
smooth, inclined boundary immersed in water or oil were evaluated from 
the relation 

„ irD 2 2 „ . 

C D p f— V s = F G Sln “ 

where v g is the measured translational velocity of the sphere of diameter D. 
Their results are shown in Figure 81. Line A in this diagram is the con- 
ventional curve for °c * for spherical particles moving in an un- 
bounded fluid (Batchelor, 1967). Eagleson et al (1958) later extended these 
results experimentally to account for the presence of a rough bottom. 

Figure 82 shows that for constant values of particle diameter D, decreasing 
roughness size k (i. e., increasing ratio D / k) causes a reduction in the drag 
coefficient determined for a smooth bottom from Figure 81. In their studies 
of incipient motion of bed particles due to surface waves, Eagleson and 
Dean (1959) chose = 19. 2/R^ with D/k =* 1 and R^ = u^ D/y. On the 
basis of a large number of experiments in which was computed from an 
equation similar to Equation (5-4) for various measured values of u^, D^, 

<j>, a, and s /s,, Eagleson and Dean (1959) concluded that there was no 

S 1 

qualitative disagreement between the values of C D obtained in this way with 
those values obtained from Figures 81 and 82. It is not obvious, however, 
that their results show conclusively that the two methods of computing 
are equivalent for incipient motion of the particles. On the other hand, 
their results are quite convincing that the two methods are in agreement 
for discrete particles that are already in motion. 

In view of the apparent importance of the local particle Reynolds 

number R^ on the value of C^, it is useful to establish some upper bound 

on R for the oceanic problem. Sediments on the outer shelf have median 
u -2 

diameters on the order of 10 cm. Lacking specific measurements of 
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FIG. 81. Resistance coefficients C for spheres (after Eagleson and Dean, 1959). 



176 




D/ k 

FIG. 82. Effect of bed-particle geometry on the resistance 

coefficients for spheres (after Eagleson and Dean, 1959). 
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internal-wave velocities in the bottom boundary layer, a value of 100 cm/sec 

* 

is selected as a crude upper estimate. 

-2 2 

If v = Vr - 10 cm /sec, then =* 100. For larger values of v 
or smaller u f , R D becomes less than 100. This upper bound suggests that, 
according to Figure 81, C D « 192/R D for a hydraulically smooth bottom. 

If the additonal assumption that D/k « 1 is made, then = 192/Rp 
(Figure 82). If this value of C n is applied to Equation (5-11) with 

U f D i 2 

R~ = and v * 5 co/2, then 

D v ' 

a w 2 

%= 31 - 36 ( 5 -! 2 ) 

for incipient downslope motion on small slopes. 

It is obvious from Equation (5-12) - and not unexpected — that the 
local incipient diameter is a function of the boundary layer thickness (or 
alternatively, the viscosity y) as well as the local depth and input wave 
conditions. The value of increases linearly with decreasing depth, sug- 
gesting progressive sorting of sediment. Equation (5-12) predicts that sed- 
iment particles having diameters smaller than will initially move down- 
slope so that there should be a trend toward smaller mean diameters in the 
seaward direction. Since the model does not describe the subsequent history 
of bottom particles once they are initially moved from their roughness hol- 
lows, it is not possible from this analysis to trace their net motion after 
initial movement. However, one might speculate that once the particles 
are set into motion they might be advected to different locations as bed load 



*A few values of instantaneous near-bottom velocities on the continental 
margins are available. For example, Emery and Ross (1968) reported a 
maximum value of 70 cm/sec on the continental slope south of Martha's 
Vineyard; the flow direction was parallel to the bottom contours. 
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if there is net motion of the surrounding fluid. The net motion of the fluid 
need not depend on mixing induced by internal waves, but can be caused by 
other currents. It is emphasized that this analysis pertains to initiation of 
motion of noncohesive bed particles for the conditions assumed; subsequent 
motion as established bed load or suspended load is only speculative. 

The inverse relationship between the ratios D./s and h/h Q can be 
conveniently illustrated for different choices of a Q , N, and u>; Figure 8? 
is an example of this relationship. Each straight line in the diagram 
represents a particular input wave amplitude a Q for internal waves that 
are assumed to be normally incident to the slope. The larger values of a Q 
are permitted only for h Q large in the sense a 0 A c ) << L (The linear in- 
ternal wave theory used in this analysis is valid only for small-amplitude 
waves.) Finite-amplitude effects have not been considered explicitly. Con- 
sequently, the results are valid upslope to a point where the Stokes param- 
eter e becomes order one; the exact value of e for this limiting condition 
is not precisely defined. A discussion of e for the small-slope approxi- 
mation is given in Chapter 4. The experiments show that for y < c, the 
high-frequency waves induce a zone of breaking for values of e > 0. 5. A 
value of e>0.5 might be used to establish a limiting depth ratio h/h Q for 
this model; i. e., at depth ratios shallower than h/h Q the linear analysis is 
invalid. If a Q = mh Q where m << 1 and positive, then from Equation (4-14) 



mode n = 1. 

For example, for relatively high-frequency, first-mode internal waves, such 
that c N/2 and a Q = 0.01 h Q (m = 0.01), it follows from Equation (5-13) 
that 




(5-13) 



where 
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FIG. 83. Ratio of incipient diameter D. to boundary-layer 

thickness S for various positions on the slope, given 
by h/h , and for selected amplitudes a Q of the 
incident internal waves. 
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If e =0. 1, h/h Q sa 0.4; if e = 0. 5, h/h Q - 0- 2. Figure 83, shows the 
limiting depth ratios derived in these examples. 

Before discussing the results of this model in terms of oceanic 
conditions it is instructive to show that the form of Equation (5-12) is in 
some sense compatible with the familiar Shields parameter S, 



duce initial motion of noncohesive bed particles in a steady flow. Shields 
(1936) presented experimental results relating this parameter to a partic- 
ular Reynolds number. Others (for example, Vanoni, 1964) provided sub- 
sequent experimental results that established a single curve relating this 
Shields parameter S to the particle Reynolds number R*^, where 



The diagram in Figure 84 shows the classical Shields curve with some 
later experimental results plotted on it; the Shields parameter is the ordi- 
nate, and the particle Reynolds number is the abscissa. Two basic prob- 
lems arise in the application of this curve to the present problem: (1) 
the experimental measurements represented by the curve were made for 
boundary layers that were fully turbulent; and (2) in most cases the mean 
flow conditions were steady (channel flow). White (1940) derived a similar 
semi-empirical result for the critical stress at the bed for steady flow 
condit ions and a laminar boundary layer: 



T 



= (s s /s f - 1) Pf g D s 



(5-14) 



where D i s equivalent to D. and t is the shear stress necessary to pro- 
S 1 c 




S = 0. 18 tan <£ 



(5-15) 



where S is given by Equation (5-14). According to White, tan <j> is about 
unity which gives S <« 0. 18, as compared to about 0.03 to 0.06 for S ob- 
tained by Shields (Figure 84). White found, however, that the value of S 
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for turbulent flow in the boundary layer was about one-half that for laminar 
conditions. The remaining discrepancies in the results of White for turbu- 
lent flow with those of Shields were attributed to the more intense turbulence 
in the flows with fully developed velocity profiles which Shields used (Vanoni, 
1964). The difference between the two results according to Vanoni was that 
for both turbulent and laminar flow conditions in the boundary layer, White 
found that r was apparently independent of R*. By contrast, Figure 84 

C D 2 

shows that r is dependent on R* to values of R* near 10 where the 
c D D 

Shields parameter approaches a constant value of about 0.06. 

The mean conditions of steady flow that are implicit in Figure 84 
are not unreasonable for low-frequency internal wave phenomena, since 
the advective accelerations are small, as was assumed in the earlier model. 
The compatibility of the earlier model with White's equation for laminar 
flow conditions (Equation (5-15)) can be demonstrated simply. If the bed 
shear stress is written as 



r = ju 9U /dr] 

and if au/a^ is approximated by a linear fit: 



au/ar? ~ Au/s » aca/c6 

1 /2 

where Au « aco/c, then for 5 = (2 v/co) ' (i. e., small slopes) and \x = v 

the bed shear stress can be rewritten approximately 



T 



V p f 



aco/c6)45/2) aco 




(5-16) 



On small slopes it was shown earlier that 

a = a h /h (correct to first order in e ). 

o o' 

Combining Equations (5-16) and (5-17) it follows that 



a o w 

2c (h/h Q ) 



(5-17) 



T 



5 



(5-18) 
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Applying this expression for bed shear stress into White's result for 
incipient motion in a laminar boundary layer (Equation (5-15)) an expres- 
sion is obtained for D./6: 

2 

a o u 

D / 6 " K gc (h/h Q ) 

(5-19) 

K = 2.78 

(p g - P f ) tan <l> 

The functional form of this result is equivalent to Equation (5-12). A closer 
comparison between this equation and Equation (5-12) is not justified because: 
(1) Equation (5-17) does not include a Reynolds number dependence; (2) 
Equation (5-19) assumes a linear velocity profile in the boundary layer; and 
(3) Equation (5-19) does not specify the point of application of the effective 
velocity acting on the bed grains. 

A brief summary of the preceding analysis might be helpful here. 

The foregoing model deals with conditions of incipient motion of noncohesive 
bed particles induced by internal gravity waves shoaling over a small linear 
slope. The static equilibrium of these particles on the bed was defined in 
terms of an incipient diameter D^. The analysis predicts that the net in- 
stantaneous hydrodynamic and gravitational moments acting on the particle 
about the contact points with other bed particles will cause those particles 
with diameters less than D. to move downslope initially. Equation (5-12) 
relates D. to the internal wave parameters (a Q , oo, c), a local depth ratio 
(h/h Q ), and the boundary layer thickness ( 6 ). It is also shown that Equa- 
tion (5-12) is compatible with the results of Shields (1936) and White (1940). 

MODEL ANALYSIS AND SPECULATION 

As defined above, the criterion in Equation (5-12) essentially 
defines the condition of neutral equilibrium set up by instantaneous hydro- 
dynamic forces of internal waves and that of gravity acting on discrete bed 
particles that are resting on a rough bed. A different dynamic balance is 
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needed to define a condition of equilibrium for particles that have already 
been placed in motion. This balance must depend on the net transport cur- 
rents influencing the particle motion, including any net currents that might 
be induced by the internal waves, and the downslope pull of gravity. How- 
ever, as was indicated in Chapter 4, steady transport of water particles 
induced by internal waves, near a sloping boundary in a stably stratified 
fluid, other than that caused by diffusion, theoretically cannot occur in a 
plane normal to the depth contours unless other processes of mixing are 
active. Consequently, it is not possible to predict a type of oscillating 
equilibrium for onshore-offshore motion of sedimentary particles for which 
D< D^. In the case of surface waves, Ippen and Eagleson (1955) and Eagle- 
son and Dean (1959) have shown that a balance between the forces on a par- 
ticle due to the mass transport current in the bottom boundary layer, directed 
upslope, and the downslope component of gravity creates a situation where 
sediment particles of increasing size in the onshore direction are in oscil- 
lating equilibrium. The diameter of the local equilibrium size is often 
referred to as the ’’null” or "equilibrium” diameter D^. The location of 
this type of equilibrium for a particular sediment size is called the ’’null” 
point. Miller and Zeigler (1964) and Eagleson and Johnson (1966) discuss 
sediment sorting by surface waves and give some comparisons with field 
data. The only available direct evidence for net fluid transport induced by 
internal waves on a slope are the laboratory results that were presented in 
Chapter 4. 

In view of the absence of other evidence for the existence of mass 
transport by shoaling internal waves, only the effects of the incipient condi- 
tion which was derived in the preceding section are considered here. Curves 
similar to that shown in Figure 83 can be used to determine the hypothetical 
distribution of along a linear bottom slope if 6 is known. Suppose that 
the initial size-frequency distribution of noncohesive bed particles is given 
by the solid curve in Figure 85 at each position along the slope. After the 
internal waves have acted to produce incipient motion of sediment having 
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Frequency 




FIG. 85. Hypothetical size-frequency distribution of 
sediments on an oceanic slope. 
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diameters less than D i at each depth according to Equation (5-12), the final 
shape of the size-frequency distribution at each location might eventually 
have a spike at and appear as shown by the broken line in Figure 85. 

The sediment size at which the spike occurs decreases in the offshore 
direction. 

It is unreasonable to require this simple analysis, with Equation 
(5-12) as its tidy result, to explain the complicated sedimentary pattern 
found on the continental margins. Instead, it is instructive to postulate 
various simple models to determine at what locations and for which hydro- 
dynamic conditions, if any, the predicted values of according to Equa- 
tion (5-12) might exceed the observed natural sediment diameters. 

Models are constructed for three provinces defined by the idealized 
bottom geometry and the simplified distribution of Brunt- Vaisala frequency 
shown in Figure 78. Each layer (i. e., the three fluid layers intersecting 
the shelf A, slope B, and rise C) with constant N is assumed to support 
small-amplitude, single-frequency internal waves of the first mode that 
propagate shoreward at approximately normal incidence to the slopes. 
Furthermore, only examples that satisfy the condition y << c were con- 
sidered, except for the hypothetical instance of waves with semi-diurnal 
frequency over the continental shelf and continental rise. In these cases, 
it is shown that the ratio y /c is on the order of one-tenth. Values of the 
pertinent hydrodynamic parameters are listed on each diagram. In addi- 
tion, since the depth range over which the analysis is valid is limited to 
small values of e (as discussed in the initial section of this chapter), the 
depth ratios h/h^at which e = 0. 1 and e = 0.5 are indicated on each diagram. 
It should be recalled that the experiments for y < c showed breaking of in- 
ternal waves at e > 0. 5. On this basis, it was suggested that the linear 
solutions might be invalid if local values of e exceed 0. 5. The sources 
for the sediment data are listed in each diagram; the approximate locations 
of the sampling areas for each source are shown in Figure 76. 
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Figure 86 treats a continental-shelf model with relatively high- 

frequency internal waves (co « N/2, N = 6 cph); the pertinent assumptions 

in the model are specified in the figure. The ordinate is mean diameter 

D m (for actual sediment data) or incipient diameter D. (computed from 

Equation (5-12)); the horizontal axis is the depth ratio h/h Q . As mentioned 

above, values of h/h Q computed from Equation (5-13) for two selected values 

of e and for the particular wave conditions (a , to, c) of the model are 

marked long the horizontal axis. The variations of D. with h/h are shown 

l ' o 

as theoretical stright-line curves for selected values of H: 



H = 



a o^ 



1/2 



A v 



o k 



1/2 



where a^ is input wave amplitude in meters, 

2 

v is viscosity in cm /sec, 

A = 1 m, 
o 

-2 2 

i/j^ = 10 cm /sec. 

cm^ • sec' 1 /^.) 



For example, the curve in Figure 86 labeled H = 5 represents 
1/2 2 - 1/2 

a v =50 cm • sec ' , or if v = v,., then a = 5 meters. This kind 

O 1C O 

of log-log plot simplifies presentation of the analytical results; the theoreti- 
cal curves are straight lines (Equation (5-12)), and can readily be obtained 

for various combinations of a and v (or 6 ). The field data of mean sedi- 

o 

ment diameters in Figure 86 show a general trend toward decreasing sizes 
in an offshore direction. This trend reverses for h/h Q >0.4 in Stetson's 
data and for h/h Q > 0. 5 in Hollister's data. 

Interpretation of the relationship between the actual mean values D 
and the theoretical distributions of D. was described earlier and will be 

l 

summarized here. D. represents the maximum sediment diameter that is in 
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ou 2 - f 2 

N 2 -oj 2 



-i 1/2 



Sediment data 
O Hollister 
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FIG. 86. Continental shelf model; relatively high frequency 
waves . 
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neutral equilibrium at slope positions given by h/h Q for hypothetical internal- 
wave types and idealized bottom conditions. The analytical model predicts 
initiation of motion of bed particles in a downslope direction if n < Ch- 
in Figure 86 incipient motion of measured sizes D is predicted for 
H > 50 (for example, a Q = 5m, v = 1 cni /sec) and h/h Q > 0.48 (for the 
cut-off condition imposed by e = 0. 5). If H = 5 (a = 5m, v = v, ) and if 

O K 

e = 0. 5 is the limiting value for the linear analysis, then all values of D 

m 

from Stetson's data for 0. 25 < h/h Q < 0.45 could be placed into motion by 
the waves. It might be noted here that the nonlinear effects which are sig- 
nificant farther upslope (where e is order one, say) could make the case for 
possible sediment movement even stronger. For instance, the experiments 
for y < c showed that the measured velocities in the boundary layer near 
the zone of breaking were significantly larger than corresponding measure- 
ments taken farther downslope. Consequently, one might expect that the 
presence of nonlinearities in the velocity field might increase the values of 
D i for positions well upslope. The consideration here of a strictly linear 
analysis is in this sense a "worst” case: if initiation of motion of measured 
diameters is predictable from the linear model, then the inclusion of non- 
linear effects strengthens the argument. We might also recall that values of 
v > 1 /^ are included mainly on the supposition that the oceanic boundary-layer 
dynamics might involve an eddy viscosity. 

Figure 87 is another shelf model similar to Figure 86 except that 
the internal-wave frequency is assumed to be semi-diurnal. In this model, 
the requirement y/c <<1 is necessarily relaxed to y/c<l (y/c == 0.01). 

The case for incipient motion of the measured sediment sizes is not as con- 
vincing as that given in Figure 86. The theoretical curves indicate that 
is less than D. for relatively small-amplitude waves (a Q < 5m) only if 
v > 1 cm^/sec (see H = 50, for example). The values of D i along theoreti- 
cal curves H = 1 and H = 5 are below the mean diameters; this suggests that 
waves of 5-meter amplitudes with v = will not move the measured mean 
sizes. 
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f = 0.062 cph 
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FIG. 87. Continental shelf model; wave of tidal frequency 
(semi-diurnal) . 
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The continental slope model in Figure 88 assumes relatively high- 
frequency waves (c o = N/2; N = 3 cph) because critical conditions (y » c) 
are approximately achieved along the slope for waves of lower frequencies, 
such as the frequency of the semi-diurnal tide. The data shown in this 
figure were obtained from various sources; the lines of sample locations are 
shown in Figure 76. Initial movement is predicted for sediment along the 
lower extremity of the slope (h/h Q > 0. 8) for H = 10. For H = 50 (a Q = 5m, 
v = 1 cm^/sec), initiation of sediment motion is possible for each of the 
measured values along the slope. 

A smaller number of data points are shown in the continental rise 
model in Figure 89. The available values of mean diameter tend to be in 
the range of fine to very fine silt; the ordinate (D^ or Dj has been scaled 
accordingly. It is worth recalling the earlier discussion of sedimentary 
types on the continental rise. Several previous workers (Heezen, et al 
1966, for example) noted that clean quartzose sands and coarse silts inter- 
bedded with the clays can be locally dominant. These larger sizes are not 
represented in the data plotted in Figure 89. It is obvious from Pratt's 
data, however, that the mean diameters found on the lower parts of the rise 
are generally larger than farther upslope. The computed incipient diameters 
exceed the mean diameters for H = 5 (a = 5m, v - v ,) on the upper rise. 
Initial movement of the larger measured sizes on the lower rise requires 
larger wave amplitudes (or larger viscosity), such that H > 50. The effects 
of cohesive forces acting on the smaller sizes (very fine silts and clays) 
have not been considered; these forces might inhibit initial movement and 
thereby require larger values of H for incipient motion. 

The results of this section can be summarized briefly. It has been 
shown that, within the constraints of the models, the predicted values of 
incipient diameters D. can exceed the values of mean diameters D mea- 
sured at various positions on the continental margin southeast of New England. 
The condition for incipient motion of bed particles induced by internal waves 
shoaling over a small slope whose average gradient is linear is given by 
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Ao= 1 ™ (average of 3 

K- 10 2 cm 2 /sec samples) 



FIG. 88. Continental slope model; relatively high-frequency 
waves . 
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FIG. 89. Continental rise model; waves of tidal frequency 
(semi-diurnal) . 
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Equation (5-12). It is obvious from Equation (5-12) that larger values of 
wave amplitude and viscosity produce larger incipient diameters at specified 
values of h/h Q . Incipient motion of the measured sediment is predicted for 
specific instances in the models considered above: 

(1) on the outer continental shelf (h > 50 m, say) with N = 6 cph, 
at = 3 cph, and H > 5; 

(2) on the outer continental shelf with N = 6 cph, semi-diurnal 
frequencies, and H > 50; 

(3) on the lower continental slope (h/h Q > 0. 8), for N = 3 cph, 

cu = 1. 5 cph and H = 10, or on the middle sections of the continental slope 
with H = 50; 

(4) on the upper and middle sections of the continental rise for 
N = 1 cph, co = 0. 083 cph, and H > 5, or on the lower rise with H > 50. 

The internal wave amplitudes suggested in the discussion of the various 
models are representative in that they fall within the range measured in 
the ocean. The various proposed values of viscosity v (or boundary layer 
thicknesses s) are speculative. 

The available data for specific sediment distributions on the conti- 
nental margins do not correspond in a simple way to the predicted distri- 
butions; the natural setting is too complex to fit these idealized models. 
However, what is not answerable at this time is whether or not the internal 
waves do in fact initiate motion of bottom sediment: for certain conditions, 
the models predict that they can. It must be noted that an additional attri- 
bute of the sediments on the slope and rise that is probably important in 
increasing the local forces necessary to dislodge the grains is the existence 
of cohesive forces acting between particles of clay and very fine silt sizes. 
The incipient-diameter criterion in Equation (5-12) assumed noncohesive 
bed particles. Other factors, such as bed roughness scales greater than 
granular, stirring by bottom organisms, bed fluidization, grain angularity, 
and packing of bed grains, are probably significant locally in the transport 
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of bed particles. Their effects are difficult to treat analytically; future 
field studies might establish their actual importance. 

In the light of the experimental results presented above, the effects 
of two other phenomena induced by internal waves must be considered. 
Breaking internal waves of high frequency were observed to produce turbu- 
lent conditions and increased instantaneous velocities near the bottom. The 
breaking often generated filimentous horizontal layers which penetrated the 
adjacent fluid interior. Secondly, experimental runs with y> c showed a 
vortex instability along the slope; the mixing in the vortices caused a ver- 
tical structure of alternating horizontal streamers and layers within which 
the net fluid transport was away from and toward the slope, respectively. 

If these processes occur in nature (there are no definitive field measure- 
ments of either process) speculation about their possible impact for sedi- 
ment processes on the continental margins raises some interesting questions. 

Wherever the breaking might occur on the slope, sediment that is 
well-sorted and coarse relative to nearby material probably will be associ- 
ated with this zone. The finer materials might be winnowed out and trans- 
ported from the area by currents. Hogg and Wunsch (1970) have shown that 
for sufficiently large values of eddy viscosity and slightly oblique incidence, 
in theory the breaking internal waves can generate "longslope” currents. 
Hence, a littoral drift of fine materials similar to that observed for surface 
waves and supplementing the effects of the contour currents (Heezen, et al 
1966) might be possible. If the position of breaking along the slope is a 
function of time, i. e., if the stratified layer containing the waves shifts its 
vertical position over a time scale of days or longer, as seasonal or local 



The laboratory measurements suggest that the breaking occurs well up- 
slope (for y < c). In fact, when the stratified section was beneath a less 
strongly stratified or homogeneous layer, the breaking occurred near or at 
the base of this upper layer. 
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thermoclines might, then a time- sequenced winnowing out of fine materials 
might occur along the bottom over a wide offshore traverse. Deposition of 
recent sediments entering the area from the continental sources might be 
prevented, and local size-frequency distributions might be skewed toward 
the coarser sizes. A dynamic process like this might help to explain the 
observed scarcity of recent sediments on the outer shelf. 

Finally, speculation is ventured for the effects of breaking internal 
waves at the edge of the continental shelf. Suppose that relatively large- 
amplitude, high-frequency internal waves were generated at infrequent 
intervals (for example, by the passage of severe surface storms) in the 
seasonal thermocline when the upper extent of the thermocline was at the 
approximate depth of the shelf edge. Sudden substantial increases in bot- 
tom turbidity due to the suspension of bed materials by breaking of these 
large-amplitude waves might increase the density of the water along the 
shelf-slope transition to values high enough to generate density currents. 
This mechanism might provide a means of continuing erosion of the slope 
face at irregular times and places by relatively weak small-scale turbidity 
currents. 

The vortex instabilities and streamers that were created in the lab- 
oratory, generally for y > c, might also provide a means of increasing 
suspension of bottom materials, by increased turbulence within the eddies 
near the bed, and by transporting fine suspended materials away from the 
slope. The net transport in the streamers might move the finer sediment 
(clay and very fine silt) away from the continental slope to deeper sites of 
deposition on the rise or abyssal plain. 

None of these mechanisms has been observed in the ocean. The 
case for increased mixing along the oceanic slopes is not new; Munk (1966) 
and Wunsch (1970) have called attention to this process. It is hoped that 
these experiments have illuminated some of the hydrodynamic features of 
shoaling internal waves, and strengthened the case for internal waves as 
geological agents on the oceanic slopes. 
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APPENDIX A 

The use of conductivity cells in A. C. circuits to detect fluctuations 
of conductivity in electrolyte solutions has been briefly summarized by 
Gibson and Schwarz (1963). In this section we describe only the particular 
sensors and circuitry that were used during the set of experiments discussed 
in the text. 

DESCRIPTION 

The sensor is a single- electrode, platinum conductivity probe whose 
construction is shown diagrammatically in Figure 4. The rounded platinum 
tip and the L shaped glass housing minimized flow disturbances induced by 
the probes. The probe assembly was oriented so that the horizontal portion 
of the glass housing containing the probe tip was approximately parallel to 
the crests of the internal waves (i. e., perpendicular to the side walls of the 
tank). This prevented the wake of the vertical shaft from washing by the 
electrode and introducing spurious values. The electrical path was between 
this electrode and the solution ground (a large copper strip placed along the 
tank bottom beneath the slope). The platinum tip was normally platinized 
by a method similar to that described by Wolbarscht et al (I960): (1) the 
probe was immersed in a 1 percent solution of Chloroplatinic acid to which 
lead acetate (80 percent by weight) had been added; (2) current from a 10 v 
D.C. source in series with a 10,000 ohm resistor was passed between the 
probe (cathode) and a platinum wire (anode) either for about 15 seconds or 
until a thin stream of bubbles was emitted from the probe tip. The result 
was a thin grayish-black coating on the platinum probe. The effects of 
this coating on probe performance are discussed in Ives and Janz (1961); 
essentially, the platinization minimizes drift in the output of the sensor. 

The probe is the active arm in an A. C. Wheatstone Bridge (Figure 
Al). Potential imbalances across the bridge caused by fluctuations of con- 
ductivity at the probe are amplified, detected (diode rectification), and 
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recorded by the A-D conversion described in Chapter 2. The high excitation 
frequency (60 kc) and low excitation voltage (0. 5 v rms) inhibited polariza- 
tion of the electrolyte at the probe tip. This was essential in order to mini- 
mize signal drift and capacitance problems. 

PERFORMANCE 

Four complete circuits like that shown in Figure AX were used. It 
was found that the levels of signal output changed by as much as 15 percent 
when probes were interchanged, owing to nonuniformities and irregular 
platinization of the probe tip; this usually necessitated a calibration check 
before each experiment. Normally this calibration consisted of recording 
output voltage as a function of probe depth before each run. This provided 
an updated sensitivity value that could be checked with previous values. 

If it was observed that sensitivity changed markedly (> 10 percent) between 
runs, the probe tip was checked for fouling and cleaned as necessary. The 
cleaning process is like that described by Ives and Johns (1961). Prior to 
initial platinization, or when the platinized coating was cleaned off with a 
soft tissue, the following procedure was used: 

(1) wash in warm concentrated HNO^ for 30 seconds, 

(2) wash in distilled water, 

(3) polish platinum tip with glass rod, 

(4) treat polished electrode with warm 50 percent solution 

of Aqua Regia (4 Volumes f^O, 3 Volumes HCC , 1 Volume 
HNOg) for 15 seconds, 

(5) wash in distilled water, 

(6) wash in warm HNO^ for 30 seconds, 

(7) wash in distilled water, 

(8) dry in air, 

(9) place in distilled water when not in use. 

Normally, the output signal was acceptable up to three hours of continuous 
use; however, shortly thereafter degradation of the signal to noise ratio 
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and erratic output values necessitated cleaning of the probe tip (steps (5) 
through (7) above). 



The following list summarizes the average performance of the 
probes determined from experimental data and calibration tests. The 
length units are presented for convenience and are based on: 



P = P Q (1 -0z) 

where 



/ 3 = 0.001 g/cc/cm 



(1) 

( 2 ) 

(3) 



(4) 



(5) 



sensitivity: 

resolution: 

accuracy: 



linearity: 



drift: 



0.3 volt/0.001 g/cc (0.3 v/cm) 

0. 5 X 10' 4 g/cc (0. 5 mm) 

5 percent of maximum vertical displacement 
of probe relative to stratified fluid at rest 
for displacements < 0. 5 cm. 

2 percent of maximum vertical displacement 
of probe relative to stratified fluid at rest 
for displacements < 0. 5 cm. 
variable, depending on quality of platiniza- 
tion. (Typically 20 millivolts/hour.) 



It might be useful to present functional relationships between some 
of the physical variables for sodium chloride solutions that were assembled 
during this study. Several sources were used to obtain the relationships 
(L.L. Higgins, 1962; International Critical Tables, 1929; American Opti- 
cal Company, 1969). 

(1) n = 1.3330 + 0.244 (p-1); T = 20°C 
(accurate to approximately 1 percent) 

(2) S = 1.39 (p-1) x 10 3 ; T = 20°C 
(accurate to approximately 1 percent) 

(3) a = 1.902 (p-1); T = 20°C 
(accurate to approximately 5 percent) 
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where 

n is index of refraction, 

3 

p is density on g/cm 
S is salinity in °/oo, 
a is electrical conductivity in (ohm- cm) 
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TABLE B1 



Group 


N 

sec -1 


N 

cps 


1 


0.776 


0.123 


2 


0.990 


0.157 


3 


0.970 


0.154 


4 


0.988 


0.156 


5 


0.945 


0.151 


6 


0.952 


0.152 


7 


0.972 


0.154 



w c -I 
sec 


w c 


T c 


cps 


sec 


0. 199 


0.038 


31.51 


0.248 


0.039 


25.35 


0.456 


0.073 


13.80 


0.494 


0.079 


12.72 


0.473 


0.075 


13.29 


0.235 


0.038 


26.69 


0.497 


0.079 


12.7 



“0 

cm 


u 

deg 


l 

(tan a) 


31.55 


14.9 


0.266 


31.60 


14.5 


0.259 


31.80 


28.0 


0.532 


31.65 


30.0 


0.577 


33.35 


30.0 


0.577 


33.45 


14.3 


0.255 


33.5 


30.0 


0.577 
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TABLE B 2 



Ref . No . 


CO 


T 


k° 


c 


q 




cps 


sec 


cm ' 1 






1-2 


0 . 129 * 


7.76 


0.142 


1.426 


17.105 


2-3 


0.129 


7.76 


0.150 


1.522 


8.614 


3-4 


0.109 


9.14 


0.095 


0.962 


4.532 


4-2 


0 . 102 * 


9.85 


0.084 


0.843 


9.913 


5-5 


0 . 102 * 


9.85 


0.083 


0.883 


4.020 


6-3 


0 . 098 * 


10.24 


0.081 


0.822 


4.075 


7-4 


0.098 


10.24 


0.079 


0.793 


3.627 


8-5 


0 . 098 * 


10.24 


0.080 


0.853 


3.818 


9-6 


0.098 


10.24 


0.080 


0.848 


10.124 


10-1 


0 . 094 * 


10.67 


0.116 


1.168 


13.546 


11-2 


0.094 


10.67 


0.074 


0.743 


8.625 


12-3 


0 . 094 * 


10.67 


0.075 


0.764 


3.656 


13-1 


0 . 094 * 


10.67 


0.116 


1.168 


13.546 


14-4 


0 . 090 * 


11.13 


0.070 


0.697 


2.897 


15-4 


0.086 


11.64 


0.065 


0.653 


2.506 


16-5 


0.086 


11.64 


0.046 


0.486 


— 


17-4 


0 . 082 * 


12.19 


0.061 


0.612 


2.085 


18-5 


0 . 082 * 


12.19 


0.061 


0.650 


2.227 


19-4 


0 . 078 * 


12.80 


0.057 


0.573 


1.537 


20-4 


0 . 072 * 


13.86 


0.051 


0.515 


-- 


21-4 


0 . 066 * 


15.06 


0.046 


0.462 


— 


22-5 


0 . 066 * 


15.06 


0.046 


0.488 


-- 


23-6 


0.062 


16.00 


0.043 


0.453 


4.922 


24-2 


0 . 062 * 


16.00 


0.042 


0.420 


4.366 


25-1 


0 . 055 * 


18.28 


0.051 


0.517 


5.518 


26-4 


0 . 055 * 


18.28 


0.037 


0.373 




27-5 


0.055 


18.28 


0.037 


0.393 


-- 


28-3 


0.051 


19.69 


0.035 


0.350 


-- 


29-4 


0 . 051 * 


19.69 


0.034 


0.343 


-- 


30-3 


0.049 


20.51 


0.033 


0.335 


-- 


31-6 


0 . 047 * 


21.33 


0.031 


0.325 


2.971 


32-6 


0 . 043 * 


23.26 


0.028 


0.296 


2.407 


33-4 


0 . 039 * 


25.60 


0.025 


0.256 


-- 
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